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1. Introduction 

The computation of amplitudes in QCD is a problem whose complexity grows factorially 
with the number of particles even at tree level. Such complexity stems from the extremely 
large number of Feynman diagrams that contribute to many-particle processes and which, 
apart from the inherent complication of the Lorentz structure, induce a proliferation of 
mutually- independent colour factors, that give rise to an involved colour algebra. Lorentz 
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and colour structures can be separated by expressing scattering amplitudes as sums of prod- 
ucts of dual (or colour-ordered) amplitudes times colour factors (see e.g. ref. [1] and refer- 
ences therein) . The problem of the efficient computation of dual amplitudes has attracted 
considerable attention, and nowadays several solutions exist (Berends-Giele recursion re- 
lations [2], CSW relations [3], BCF recursion relations [4]). In order to predict observable 
cross sections, however, the colour algebra must be performed. Its factorial growth has 
been bypassed by working in the colour-configuration space [5], and by re-expressing it 
as an ordinary integral [6] ; these approaches avoid the use of Feynman diagrams [7] , and 
are essentially equivalent to colour-dressed recursion relations [8]; by these means, the 
Feynman-diagram factorial complexity is reduced to an exponential one. 

The problem posed by the factorially-growing complexity of the colour algebra was 
simply irrelevant for the calculation of observables at the next-to-leading order (NLO) and 
beyond, because several other issues limited anyhow the applicability of NLO techniques 
to small-multiplicity processes. This is not the case any longer, thanks to the progress 
recently achieved in the automation of the two essential ingredients in an NLO computa- 
tion: the subtraction of real-emission singularities [9, 10, 11, 12] (where use has been made 
of the universal subtraction formalisms known as FKS [13, 14] and dipole [15]), and the 
computation of one-loop matrix elements [16, 17, 18, 19, 20] (based on generalized unitar- 
ity [21, 22, 23] and integrand reduction [24, 25, 19] techniques). Therefore, the question 
becomes relevant of how to best organize an NLO computation, in order to be able to ex- 
ploit the solutions that work well for tree- level matrix elements, and to use new tree-level 
approaches which may become available in the future. As far as the one-loop contribution 
is concerned, the problem is conceptually analogous to that relevant to tree-level ampli- 
tudes (but technically more complicated, see e.g. ref. [26]), since it consists in finding the 
optimal representation of an amplitude in terms of dual amplitudes and colour structures. 
On the other hand, the subtraction of real-emission singularities is complicated by the fact 
that it must be performed at the level of amplitudes squared, while it is the amplitudes 
(not squared) that display the simplest factorization properties in the soft and collinear 
regions (see e.g. ref [27]). As a result, in the subtraction terms defined in the FKS and 
dipole formalisms, the colour structure is not factorized, since it appears in both the uni- 
versal kernels and the short-distance, process-dependent reduced matrix elements. Such 
a convolution is avoided in the antenna subtraction method [28, 27, 29], which is based 
on the use of the squares of, or the interferences between, dual amplitudes as elementary 
quantities whose singularities have to be subtracted. 

The aim of this paper is that of formulating the FKS subtraction in terms of colour- 
dressed and colour-ordered amplitudes. This will allow one to use, with only a few triv- 
ial modifications, the techniques developed at tree level to deal with the problem of the 
factorially-growing complexity. In order to be more precise, let me briefly digress and 
introduce the tenets of the FKS method. 

The basic idea of the FKS subtraction formalism is that of treating in an independent 
manner the singularities present in a multi-parton matrix element squared. This is achieved 
by damping all singularities except one soft and one collinear, which are pre-determined, 
and by repeating this procedure for all singularities in turn. One introduces a set of 
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functions Sij such that: 



J>, = 1. (1.1) 



The sum in eq. (1.1) can be thought of as extending to all strongly-interacting pairs of 
particles, although in practice significant simplifications are possible (see ref. [11] for an 
exhaustive discussion). For a given pair, the function Sij is equal to zero in all soft and 
collinear limits, except those associated with ki — > 0, and with k{ || kj. Given a real-emission 
matrix element squared M, one exploits eq. (1.1) by writing 

M = Y^ M ij > Mij = S i:j M . (1.2) 

The matrix elements Mij are independent from each other and, thanks to the properties 
of the Sij functions, have one soft and one collinear singularity at most (depending on the 
identities of particles i and j). These singularities are subtracted as follows: be Ei and 
Qij the energy of parton i and the angle between partons i and j, in the rest frame of the 
incoming particles. The (divergent) integration of Mij over the phase-space d(f> is replaced 
by its (convergent) subtracted form: 

Miji4 , (1) + ( T -i^) + [£?(! -co S% )M,] f . (1.3) 

If one neglects, for the sake of this discussion, the factors that multiply M in the square 
brackets on the r.h.s. of eq. (1.3), one sees that the result of the two plus prescriptions is 
that of constructing the linear combination: 

M - M SOFT - M COLL + -MsC • (1-4) 

Here, the second, third, and fourth terms are the limits (in the sense of asymptotic be- 
haviour, as customary) of M when Ei — > 0, 9^ — > 0, and (Ei,0ij) — > (0,0), respectively. 

As was already mentioned, so far the FKS procedure has been formulated by assuming 
that all contributions to eq. (1.4) be colour-summed. In the present paper, M will be 
decomposed into a sum of terms, each of which has an immediate interpretation in terms 
of colours. Furthermore, this sum will have to commute with the subtraction procedure. 
In other words, each of the terms in the sum must have well-defined soft, collinear, and 
soft-collinear limits, so as for each such individual term a linear combination identical to 
that of eq. (1.4) can be defined, which is finite locally in the phase space. I shall consider 
two options: the sum over colour configurations, and the sum over colour flows. In the 
former case, the relevant matrix elements will be the colour-dressed amplitudes squared; 
in the latter case, they will be the squares of, or the interferences between, colour-ordered 
amplitudes. 

This paper is organized as follows. In sections 2 and 3 I shall carry out the programme 
described above for pure-gluon amplitudes and for generic quark-gluon amplitudes respec- 
tively. The reader who is not interested in the technical details of the derivation may 
skip these two sections, and go directly to section 4, where I summarize the results for 
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the limits of the matrix elements. In section 5 I consider the case of the contributions to 
the NLO cross section which do not originate from the subtraction procedure of eq. (1.3). 
Section 6 presents a discussion on the findings of this paper and a very brief comparison 
with antenna and dipole methods, and section 7 reports my conclusions. The conventions 
adopted for the colour matrices are given in appendix A, while appendix B contains a few 
technicalities, whose role is important for the understanding of the derivations presented 
in the main text. 

2. Gluon amplitudes 

Given that the aim of this paper is the re-organization of the subtraction procedure in 
a colour-friendly way, it is convenient to consider all particles entering a hard scattering 
as outgoing, in order to simplify the notation. The Born and real-emission processes will 
therefore be 

— >n, (2.1) 
— >n + l, (2.2) 

respectively. Furthermore, it is not restrictive to assign labels n+1 and n to the FKS parton 
and to its sister, respectively 1 . The final results which will be obtained here will be easily 
cast into the usual form adopted in the context of FKS subtraction elsewhere [13, 14, 11] 
- processes with a physical four-momentum configuration 2— > n — 2 or 2— > n + 1 — 2 will 
be obtained by crossing, while partons (and in particular the FKS parton and its sister) 
may be relabeled. In the introductory part of this section, where I shall discuss the general 
features of the scattering amplitudes, I shall consider the Born process of eq. (2.1) to be 
definite. The case of the real-emission process of eq. (2.2) can simply be obtained with the 
formal replacement n — > n + 1. 

Given an n- gluon colour configuration in SU(./V) 

{ai}?=i, <k e {l,...iV 2 -l} (2.3) 

the corresponding scattering amplitude can be written as follows, where I adopt the rep- 
resentation used in ref. [1]: 

^(oi, . . . an) = Tr(AMD . . . \ a ^A (n) ((j(l), . . . a{n)) 

ee £ A ({a*}, a) (2.4) 

where I have introduced the shorthand notation: 

A({ ai },CT) = Tx(\ a °m ...X a ^y (2.5) 

X I remind the reader that, within one given region of the FKS dynamic phase-space partition achieved 
by means of eq. (1.1), the FKS parton is defined to be the only one that can give rise to soft singularities 
(hence, i), while the pair composed of the FKS parton and its sister is the only pair that can give rise to 
collincar singularities (hence, 
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When no ambiguity is possible, the first argument of A (i.e., the set of colour indices) will 
be dropped. Here, A a are the Gell-Mann matrices, whose normalization conventions are 
given in appendix A, and P' n is the set of non-cyclic permutations of the first n integers. 
The non-cyclicity condition can be imposed by simply requiring 

a(l) = 1. (2.6) 

The quantity on the l.h.s. of eq. (2.4) is the colour-dressed amplitude, and A^ is the dual 
(or colour-ordered) amplitude. Each term in the sum on the r.h.s. of eq. (2.4) corresponds 
to a colour flow (in which the colour of gluon a(i) is connected with the anticolour of gluon 
a(i + 1)). In a technical sense, I shall identify the flow with the ordered set 

Kl),...a(n)) . (2.7) 

It will be convenient to regard each colour configuration as a vector \ai, . . . a n ); the set of 
all colour configurations can therefore be made equivalent to an ortho-normal base in a 
vector space (called colour space henceforth): 

n 

(h, ■ -AK, ...a n ) = Y[6 ai b t , (2.8) 

i=i 

^2 \a 1 ,...a n )(a 1 ,...a n \ = I. (2.9) 

K}?=i 

Equation (2.4) suggests to define the following vectors in the colour space: 

\A^( ai , . . . a n )) = A^( ai , . . . a n )\ ai , ...a n ), (2.10) 
\A^ n \a)) = A({a i },<7).AW((7)|a 1 ,...a n >, (2.11) 

which can be used to construct a vector that corresponds to the physical amplitude (i.e. 
the amplitude obtained by summing over all colour configurations or flows): 

\A™) = \A {n \ ai ,...a n )), (2.12) 

{Oi}?=l 

= l-4 (n V)>- (2.13) 

What is relevant to cross-section computations is the amplitude squared 2 : 

M {n) = (A {n) \A {n) ) . (2.14) 

Equations (2.12) and (2.13) can be used to rewrite the amplitude squared in two different 
ways: 

M (n) = £ M^( ai ,...a n ) (2.15) 
= Y, X (n V», (2-16) 

*,a>eP> 



2 In refs. [13, 14, 11], the notation is used for an amplitude squared, times the flux factor, times spin 

and colour average factors. All these factors are omitted here, being irrelevant for the present discussion. 
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where I defined 



M W (ai, ■ ■ ■ a n ) = {A&Hau . . . a n )\A (n \ ai , ■ ■ ■ O) (2.17) 
= A^( ai ,...a n )\ 2 , (2.18) 
M (n) (a',a) = (A {n) (a')\A (n) (a)). (2.19) 

In terms of scalar quantities, eq. (2.19) is: 

M^(a',a)=A (n \a , yC(a',a)A^(a), (2.20) 
where I defined the colour- flow matrix element: 

C(a',a)= £ A(KU')*A(KU) . (2.21) 

The hermiticity of the Gell-Mann matrices and eq. (A. 5) imply that C is a real, symmetric 
matrix: 

C(a,a')eR, C(a',a) = C{a,a'). (2.22) 

Note that while the amplitudes at fixed colour configurations are orthogonal (a condition 
which is formally enforced here by eq. (2.8)), this is only true for |„4,( n )(<r)) in the large- N 
limit. In fact, it is well known that 

M (n \a',a) = c a 5 a/a N n ~ 2 (N 2 - l) + O (J^j , (2.23) 

with c a a suitable real number. From the definition, one also sees that 

M (n) (a',a) =M (n) (a,a'y. (2.24) 

The quantities defined in eqs. (2.17) and (2.19), or their analogues for quark-gluon 
amplitudes to be introduced later, will be the basic building blocks for the definition of 
FKS subtraction at fixed colour configurations and flows respectively. As far as the latter 
is concerned, a final comment is in order. Colour flows are an important ingredient in 
the context of event generators, where they are used to determine, on statistical basis, the 
colour connections amongst the hard partons which initiate the showers. This determina- 
tion is driven by Ai^ n \a', a) with a' = a (since such quantities are positive-definite), while 
the matrix elements with a' ^ a are taken into account only in an averaged sense [30]. 
Therefore, as far as event generators go, a flow may actually be better defined as the pair of 
identical permutations (a, a), the pair understanding a quantity relevant to the amplitude- 
squared level. In order to generalize this idea in a way consistent with the terminology 
typically used when dealing with amplitudes, I shall call the pair 

(a', a) (2.25) 

a closed flow. This reminds one of the fact that the pair in eq. (2.25) corresponds to a set 
of colour loops; the counting of loops is a very simple way to estimate the largest possible 
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power of N which appears in A4^ n \a',a). To further this, as one sees from eq. (2.21) the 
two flows that define a closed flow play a similar, but not identical, role. When I shall need 
to distinguish them, I shall call the ordered sets 

(a'(l),...c/(n)) , Kl),...a(n)) (2.26) 

as L-flow and R-flow respectively, with the understanding that they enter the amplitudes 
{A( n \<j')\ and |„4,( n )(cr)). This naming convention stems from the simple interpretation of 
eq. (2.19) in terms of cut diagrams, where one (arbitrarily) associates bra vectors (.| with 
the left side of the cut. 

2.1 Subtraction of colour-summed matrix elements 

According to the conventions introduced at the beginning of sect. 2, the quantities that 
enter eq. (1.4) are defined as follows: 

lim M^ n+1) = Mott , (2-27) 

fc„+l->0 

lim M (n+1) = MoIl • (2.28) 

k n + l \\k n 

The counterterm A^sc +1 \ which is responsible for removing the double counting due to 
A4soft^ and AIcoll? m the soft-collinear region, is by construction: 

A4c +1) = lim -Msoft = lim M<£$ . (2.29) 

k„+i\\k n k n +i^0 

The last identity in eq. (2.29) can, and will, be used as a check of self-consistency when 
constructing A^oft* an d m'cou! ■ 

The result for the soft limit of real-emission matrix elements can be taken e.g. from 
ref. [11]: 

1 n 

M^ = \g 2 s Y. MM?; (2-30) 
k,l=l 

with 

kk-h 



kJ = Ti- !L iri— (2-31) 

Kk ' K n+1 K l ' 

M% } = -2{A<ri | Q b ( k )Q b ( l ) \ A(n) ) ■ (2-32) 

b 

The term in round brackets on the r.h.s. of eq. (2.31) is redundant, since k\ = 0, but it 
is useful when one formally manipulates eikonal factors. The operators Q b (k) describe the 
way in which the colour of gluon k is affected when a soft gluon (which here is always 
labelled by n + 1) of colour b is emitted by it. In the present context, it is useful to regard 
the result of Q b (k) as twofold: it changes the colour state of gluon k, and it creates gluon 
n + 1 with colour b. This can be written as follows: 

(a k a n+1 \Q b (k)\c k ) =6 ban+1 (Q b (k)) = 6 ban+1 (V) , (2.33) 
V) = -if bac . (2.34) 
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The definition given here implies that Q b (l)\A^) is a vector in the colour space of n + 1 
gluons. In the colour subspace of gluon k, the operator Q b (k) is hermitian. 

The collinear limit of real-emission matrix elements can also be taken from ref. [11]. 
It reads: 

M — ~ A{ p ^ )Min)+Q ^H^ M - 1 ) I (2 ' 35) 

where 

M^ + = {A^\A^), (2.36) 

with A^± the n-gluon amplitude at fixed ± helicity of the n th gluon (i.e. of the gluon that 
branches). The momentum fraction z is defined according to: 

k n = z(k n + k n+l ) when k n+l \\ k n . (2.37) 

The symbols P gg {z) and Q gg *(z) denote the Altarelli-Parisi kernel and its azimuthal coun- 
terpart, respectively, relevant to g — > gg branchings. The latter kernel has been introduced 
in ref. [13], and I refer the reader to that paper for further details. By using the identity 

Y,Q b {n)Q b {n) = C A I , (2.38) 

b 

with / being the identity operator in the colour space, I now rewrite eq. (2.35) in a slightly 
different (but completely equivalent) form, which is better suited to the manipulations I 
shall carry out in the rest of this paper: 

M - - dM^"" }■ (2 - 39) 

where 



M W = ^ (n) | Q\ n )Q h {n) \A in) ) , (2.40) 
b 

M { l l l = (A™\ ^Q b (n)Q h (n) \A?) , (2.41) 



b 

and I have introduced the quantities 



Pgg(*) = ^rPgg(z) , QgA*) = ^Qgg< z ) ■ ( 2 -42) 

^ A ^ A 

The consistency between eqs. (2.30) and (2.39), in the sense of eq. (2.29), is easy 
to prove. Let me start from computing the collinear limit of A^soft^- Since the only 
dependence on k n+ \ is in the eikonal factors, the relevant quantity is: 



lim 



1 ; (Skn + Sm)(l-S kl ), (2.43) 



1 z k n ■ k n +i 



-8- 



having used eq. (2.37). Therefore 

2 i if n ~ 1 n ~ 1 1 

am M?«> = fr^r^- E-< + E*£? • < 2 - 44 > 

As is proved in appendix B, the colour-conservation condition is 

n n—l 

Y J Q\k)\A {n) ) = ^Q 6 (fe)|^( n )) = -Q 6 (n)|^W) (2.45) 

fc=i fc=i 

for any colour index 6. Therefore, by using the definitions of the colour-linked Born's given 
in eq. (2.32), and the ortho- normality of the colour- vector basis of eq. (2.8), one obtains: 

lim a4ott } = ff L 2 , I (A^\ ^2Q b (n)Q b (n) \A™) . (2.46) 



This expression is manifestly identical to the soft limit of eq. (2.39), since 

^W^Y^. ^ 0. (2.47) 

A comment, which will apply throughout the paper, is necessary here. In general, the 
reduced kinematic configurations (that enter the Born amplitudes) that one obtains by 
taking the soft or the collinear limit of a fully-resolved configuration (that enter the real- 
emission amplitudes) need not coincide. However, as explained in ref. [11], in the context 
of FKS subtraction is possible and convenient to adopt phase-space parametrizations for 
which the two do coincide. This justifies the fact that I have used the same symbol A^ 
e.g. in eqs. (2.32) and (2.40), which would in general understand different kinematics. I 
shall use a unique notation for the soft- and collinear-induced Born amplitudes in the rest 
of the paper. Having said that, I should also like to stress that, even if the two kinematics 
were different, the results would be unchanged, since in the FKS method soft and collinear 
singularities are treated separately. The only exception would be in sect. 6, where the 
choice of kinematic configurations in the two limits becomes relevant. 

2.2 Subtraction at fixed colour configurations 

In this section, I shall consider the problem of defining the subtractions of eq. (1.4) for a 
given colour configuration. This implies that the quantity whose soft and collinear limits 
one needs to construct is 

M^ n+1 \ ai , . . . a n+1 ) = |-4 (n+1) (a!, . . . a„ +1 )| 2 . (2.48) 

The relevant limits can be obtained directly from A^soft^ an d At^JoL^ i given in eqs. (2.30) 
and (2.39) respectively. Formally, one proceeds as follows. One starts from eqs. (2.27) 
and (2.28), and then fixes the colour configuration {aj}™^ 1 on both sides of those equations 
by inserting there the projector 

\a 1 ,...a n+1 ){a 1 ,...a n+1 \. (2.49) 
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In the l.h.s.'s, such an insertion indeed results in singling out the colour-dressed amplitude 
squared of eq. (2.48): 

M (n+1) = (^(n+l^Cn+l)) 

— »• <-4 (n+1) |ai, • • • an+iXai, . . . a n+1 |^ n+1 )> = M^ n+1 \ ai , . . . a n+1 ) . (2.50) 

When performing the insertion of the projector of eq. (2.49) in the r.h.s. of eq. (2.27), by 
using eq. (2.30) one obtains: 

Mttl\a 1 ,...a n+1 ) = -glJ2 [ k > l ] M?K . . . a n+1 ) , (2.51) 

k,l=l 

where 

M%\ ai , . . . a n+1 ) = -2(A^\ Q a "+i(k)\ ai , . . . On+iXd, . . . a n+1 \Q a ^(l) \A™) , (2.52) 

having taken the Kronecker delta of eq. (2.33) into account. Using eqs. (2.12) and (2.33), 
eq. (2.52) can be easily re-expressed in terms of scalar quantities best suited to numerical 
computations: 

M% ) (a 1 ,...a n+1 )= (2.53) 
- 2^2A {n) (a 1 ,...b' k ...ai . . . a n )*Q gg (a n+1 ; a k , a^y^A^ {ai, . . . a k ...b t ...a n ), 

where I have introduced the matrices 

Q gg (a n+1 ;a k , ai ) bc = (T^) 6(lfc (T«"+% c = f^+^b f a n+1 a lC {2M) 

As the notation suggests, I stress that the colour index a n +i is not summed over on the 
r.h.s. of eq. (2.54). By construction, in SU(3) there are 8 3 Q gg matrices, of dimension 8x8. 
The computation of the collinear limit is performed along the same lines. One obtains: 

Moll 5 («!,••■ «n+i) = Tr ^—\p gg (z)M ( - n \a 1 ,...a n+1 ) 

+ Qgg*(. z ) 3? ^|^|M n |(a 1 , . . . a n+1 )) |, (2.55) 

with 

M {n \ ai ,...a n+1 ) = (A^\Q a ^(n)\a 1 ,...a n+1 )(a 1 ,...a n+1 \Q a ^(n) \A^) , (2.56) 

M n |(ai,...a n+1 ) = (A H \Q a "+ 1 (n)\a l ,...a n+l )(a 1 ,...a n+1 \Q a ^(n) \A { ?) . (2.57) 

These can be re-expressed in terms of scalar quantities, as done for their soft counterparts 
in eq. (2.53): 

M n) (ai,...a n+ i) = 

^2 „4 (n) (ai, . . . a n -i,b' n y Q gg (a n+1 ; a n , a^b^A^ (a l7 . . .a„_i,6„), (2.58) 

M^liai, . . . a n+1 ) = 

^2 A^\ai, . . .a„_i,6' l )*Q ss (a ri+ i;a„,a n )^ f , n ^l ( ( ra) (ai 1 • ■ • a n _i, &„). (2.59) 
b' n b n 
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The consistency between eqs. (2.51) and (2.55) in the sense of eq. (2.29) can be proved in 
the same way as was done at the end of sect. 2.1 for eqs. (2.30) and (2.39). In fact, apart 
from eq. (2.43) which holds independently of the treatment of the colours, the key point 
there is colour conservation as given in eq. (2.45). As is discussed in appendix B, colour 
conservation is a property that holds at fixed (n + l)-gluon colour configurations, which is 
sufficient to prove the point. Equations (2.54) and (A. 4) imply that 

£ £ Qgg( a n+ii a n , a n ) bc = C A 5 bc , (2.60) 

which is nothing but the matrix form of eq. (2.38). 
2.3 Subtraction at fixed flows 

I now turn to the definition of subtraction at given flows. The problem is formally more 
involved than that of sect. 2.3, because of the lack of the analogue of the projector in 
eq. (2.49). The basic idea is however still the same. To be definite, I consider the soft 
limit, the collinear limit being fully analogous. One starts from eq. (2.27), and expresses 
both sides as sums over (n + l)-gluon L- and R-flows: 

lim I V A^ (ri+1) (£',£) | = V lim A^ (n+1) (£', £) 

= £ Mof^ (£',£)• (2.61) 

At this point, since the representation in terms of flows is unique, the terms with the same 
£ and £' on the two sides of eq. (2.61) will be the analogues, at fixed flows, of the first two 
terms of eq. (1.4) - they will play the same roles as eqs. (2.48) and (2.51) played in the 
case of subtraction at fixed colour configurations. 

This procedure is completely trivial as far as the l.h.s. of eq. (2.61) is concerned, since 
it amounts to using the definition of colour vectors at fixed flows, given in eq. (2.19): 

^(n+i)(S',E) = (A {n+1) {V')\A {n+1) {X)) 

= A {n+1) (£')* C(£', E) A (n+1) (£) , (2.62) 

with the colour-flow matrix C defined as in eq. (2.21), with n — > n + 1 there. What is 
less trivial is the r.h.s. of eq. (2.61). The existence of such a decomposition is guaranteed 
by the fact that the amplitudes relevant to the reduced matrix element A4 SOFT live in 
the (n + l)-gluon colour space. However, their colour structures are obtained by means 
of the operators Q b that act on n-gluon colour vectors (see eq. (2.32)), and one therefore 
obtains the contributions at given £ and £' in a rather indirect way. The easiest way to 
perform the computation of A^sopt^E', £) is obviously that of using the decomposition of 
the rt-gluon colour vectors in terms of n-gluon flows, eq. (2.13), as the starting point. This 
implies that when computing A4g OFT one will end up dealing with quantities such as: 

£ [k,l] (A^(a')\ ^2Q b (k)Q\l) \A^(a)). (2.63) 

k,i b 
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As is shown in appendix B (see in particular eq. (B.10)), for a given index I and an n-gluon 
R-flow a, the vector Q b (l)\A^(a)) corresponds to the two (n + l)-gluon R-flows 

I + {a-\l))a : I-(a-\l))a, (2.64) 

with I± defined in eqs. (B.7) and (B.8). In the case of L-flows, one finds instead: 

I + (a'~\k))a' , I-(a'~ 1 (k)) a' . (2.65) 

This is equivalent to saying that eq. (2.64) defines two maps: 

(j ><7 ) i±> S + , (I, a) ^ E_, (2.66) 

relevant to R-flows. The situation is obviously identical for L-flows: 

(k, a') E' + , (k, a') ^ E'_ . (2.67) 

In order to find the representation that appears on the r.h.s. of eq. (2.61), one has therefore 
to fix E± and E' ± , and to invert the maps of eqs. (2.66) and (2.67). 

I shall explicitly carry out this procedure in sect. 2.3.1. Before turning to that, it is 
worth remarking that by construction the gluon-insertion operators Q b do not affect the 
dual amplitudes (i.e., they leave the Lorentz structure invariant). Hence, both Q b \A (jl \a)) 
and \A^ n \a)) will contain the same n-gluon dual amplitude A^ n \a). This suggests a 
procedure alternative to that implied by eq. (2.61). Namely, one may simply fix the n- 
gluon flows that appear in eq. (2.63), and therefore choose not to invert the maps of 
eqs. (2.66) and (2.67). Rather, these maps will be used to identify all M^ n+1 \T! , E) whose 
sum has a soft limit proportional to the quantity in eq. (2.63). I shall discuss this approach 
in sect. 2.3.2. 

The physical meaning of the two procedures sketched above is obvious. While in 
eq. (2.61) one fixes the flows at the level of real-emission matrix elements, in the other case 
the fixed flows are those of the Born matrix elements. Therefore, I shall refer to these two 
viewpoints as fixed real flows and fixed Born flows respectively. 

2.3.1 Fixed real flows 

As was anticipated above, I start by expressing the soft matrix element as a sum over 
n-gluon flows: 

a4ottV», (2-68) 

1 n 

-glJ2 M]<m£V,*)> (2-69) 
k,l=l 

-2{A^(a')\ Q b (k)Q b (l) \A^(a)) . (2.70) 

b 



M (n+1) 

- /vt SOPT 
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In order to compute explicitly the colour-linked Born's that appear in eq. (2.70), I introduce 
the following (n + l)-gluon colour vectors associated with an underlying n-body dynamics: 

|4x ) t(S))= Yl A({a i },E)^")(E lH4 .)|ai,...a n+ i) ) (2.71) 

where S G P^ + i, and I defined 

E a+3r =(E(l),...jH^:,...E(n + l)) . (2.72) 
By using the results of appendix B, and as already anticipated in eq. (2.64), one obtains: 
£Q b (Z) \A^(a)) = \A^ T (I + (a-\l))a)} - \A^ T (I-(a-\l)) a)). (2.73) 

b 

I can now use eq. (2.73) in eq. (2.70) since, for any two n-gluon colour vectors (v'\ and \v), 
the follow identity holds: 

J2(v'\Q\k)Q\l)\v) = ^5bc(v'\Q\k)Q c (l)\v) =Y,Sbc{w'b\wc) = ^2(w'b\wc) , (2.74) 

b b,c b,c b,c 

since 5t, c is automatically enforced by the dot product of the colour vectors, which receives 
the contribution (b\c) = 5b c from the (n + l)* ft -gluon subspace. Then: 

(</,*) = -2{(A^ T (I + (a'-\k))a')\ - <4^t(/-(^ _1 (A;)) a')\] 

{\A& T (I + (<J-\l))a)) - |^ t(j _ (<t -i (z))<7)) j . (2 . 75) 

One can now change the labellings in the sums that appear on the r.h.s. of eq. (2.69), 
J2klf(k,l) = Ylij /( (7 '(*)) <T (i))- Furthermore, one exploits eq. (B.13), which can also be 
conveniently extended: 

J_(l)<7 = J+(0)<7 , |^ T (J + (0)<7)) = \A^ T (I + (n)a)) , (2.76) 

where the last condition can be imposed because of the cyclicity of the trace. One then 
obtains: 



i 

MotxVv) = -9 2 s E <^t(/+(*)o / )i^ t (/ + o») 

- {A^ T (I + (i) a')\A^ T (I + (j - 1) a)) - {A^ T (I + (i - 1) a')\A^ T (I+(j) a)) 

+ (A%> T (I+(i - l)a')\A<£ T (I + (j - l)a)> J . (2.77) 

By relabelling j — 1 — > j in the second and fourth terms on the r.h.s. of eq. (2.77), and 
i — 1 — >■ i in the third and fourth terms, and by using the second equality in eq. (2.76), one 
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finally gets: 

- a'(i + l),a(j)\ + [a'(i + l),a(j + l)\} 

x <^ T(/+(i) a ')\ A ^ T (I + (j) a)), (2.78) 

where, as a consequence of the relabellings mentioned above, one must understand 

a(n + 1) = cr(l) . (2.79) 

At variance with eq. (2.70), eq. (2.78) contains a single (n + l)-gluon closed flow for a given 
pair. This is what renders it easy to cast the sum over a and a' (see eq. (2.68)) of the 
quantities in eq. (2.78) in the same form as the r.h.s. of eq. (2.61). One can start by using 
the following identity, which I write for R-flows to be definite (the case of L-flows being 
identical) : 

n n 

E E/(^') = E E E S(^,I + (j)a)f(a,j), (2.80) 

with 5 on the r.h.s. of eq. (2.80) the Kronecker delta. The identity above holds true since 
1= ^(S,J + (j» forgiven a, j , (2.81) 

SGP n+1 

which is a consequence of the fact that, at fixed a and j, the (n + l)-gluon flow I+{j)cr 
always exists and is unique. Now, the sum over S in eq. (2.80) will play the same role as 
that on the r.h.s. of eq. (2.61). Hence, one can exploit the 5 in eq. (2.80) to get rid of the 
sums over a and j. This implies solving the equation: 

/+(?> = £ (2.82) 
(a(l), . . . a(j),n + l,a(j + 1), . . . a(n)) 
= (S(l), . . . S(j), S(j + 1), S(j + 2), . . . E(n + 1)) (2.83) 

for a given S. It is immediate to see that this solution always exists: 

o = (2-84) 
j = E- 1 (n + l)-l, (2.85) 

and is unique. In other words, there is one-to-one correspondence between pairs composed 
of one P' n permutation and one integer, and P' n+ i permutations. Equation (2.80) then 
becomes: 

n 

EE/M= E + ■ (2-86) 

aeP^ j=i seP n+1 
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Equation (2.86) is consistent with simple counting: the n\ P' n+ \ permutations can be 
obtained by inserting the integer n + 1 in all the n possible ways (not n + 1, because of 
non-cyclicity) in each of the (n — 1)! P' n permutations: n! = n x (n — 1)!. 
These arguments, when applied to L-flows, give 

a' = E^r, (2.87) 
i = S /_1 (n + 1) - 1. (2.88) 

The results above can now be used in eqs. (2.68) and (2.78). One obtains: 

.(n+l) _ A//("+ 1 ) 



-M^t } = £ (2.89) 



s,£'ep; +1 



(n+l) 



(£',£) = - 5 |{[s'(E'" 1 (n + l)-l),S(S- 1 (n + l)-l)] (2.90) 
S'(S /_1 (n+l) - 1), S(S~ x (n + 1) + 1) 
S / (S /_1 (n + l) + l),S(S- 1 (n + l) - 1) 
+ [S'(S /_1 (n + 1) + 1), S(S- x (n + 1) + 1)] }m^ t (S', £), 
M^t(S',E) = (^t(S')I^t(S)). (2.91) 

In the evaluation of the eikonal factors that appear in eq. (2.90), one must take the ana- 
logues of eqs. (2.6) and (2.79) into account, namely 

S(l) = l, £(n + 2) = £(l), (2.92) 

and similarly for S'. Note that the second identity in eq. (2.92) is enforced by the solution of 
eq. (2.83) in the case when n+l occupies the rightmost position in I+(j)a. Equations (2.90) 
and (2.91) only depend on S and and are therefore in the form suited to be used in 
eq. (2.61). In terms of scalar quantities, eq. (2.91) reads as follows: 

M^ T (S',S) = >)(S^C(E',S)^)(S^). (2.93) 

I now address the case of the collinear limit, which I deal with as was done for the soft 
limit. Namely, in eqs. (2.39)-(2.41) I use the representation of the amplitude in terms of 
flows, eq. (2.13), and its analogue for L-flows: 

(^W| = (^ (n V)l- (2-94) 

In doing that, I explicitly express the real part that appears in the term proportional to 
Q gg * in eq. (2.39), by writing it as the sum of its argument plus the complex conjugate of 
the latter; this is necessary in order to identify a with the R-flow, and a' with the L-flow. 
I thus obtain: 

^ffi= E Mol'V>), (2-95) 
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where 

Mo + llV>) = Tr^\P 99 {z)M {n) W,a) (2.96) 

_i_ (*,\ I (k n k n +l) -jj(n) , , s [k n k n+ i] -jrjjn) , , s 

Z \ [K n K n+ i\ \Kn K n+l) 



and 

7W (n) (a',a) = <^V)I ^Q b (n)Q b (n) \A^(a)) , (2.97) 

b 

M { £(a',a) = <4V)I ^Q b (n)Q b (n) \A^\a)) . (2.98) 

( n \ 

Given that eq. (2.98) is essentially identical to eq. (2.97), I shall deal only with Ai in 
the following; the term proportional to the azimuthal kernel Q gg * (z) will be reinstated at 
the end. By comparing eq. (2.97) with eq. (2.70) one gets: 

M (n \a',a) = -^M^(a',a). (2.99) 
Hence, from eq. (2.75) one obtains: 

M {n \a',a) = (A^UlA^\n))a')\A^Ul + (^Hn))a)) 

- (A™ T {I+{o'-\n))</)\A™ T (I-{*-\n))<T)) 

- (A^(I-(</-\n))</)\J$ T (I + {*- 1 (n))*)) 

+ (A^ T {I-{a^\n))</)\A^ T (I-{a- 1 {n)) < r)). (2.100) 

The (n + l)-gluon L- and R-flows that appear on the r.h.s. of eq. (2.100) are such that: 

I + {a'~\n))u' , /+( C 7- 1 (n))c7 -> (. . . n, n + 1, . . .) , (2.101) 
I-(o'~\n))a' , I-(c- l (n))a — > (. ..n + l,n,...), (2.102) 

that is, n and n + 1 must be contiguous. This implies that all flows for which this condition 
is not satisfied correspond to matrix elements that do not diverge in the collinear limit. In 
order to proceed as was done in the case of the soft limit, I can now impose the analogues 
of eq. (2.82). By considering only R-flows to be definite, there are two possible 
shown in eqs. (2.101) and (2.102), and therefore one has either 

J + (ff _1 (n))o- = S = (. ..n,n+l,...), (2.103) 

or 

I-ia- 1 ^))^ = S = (...n + l,n,...). (2.104) 
These equations impose the following constraints on £: 

£ _1 (n) = £ _1 (n + 1) - 1 (2.105) 
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in the case of eq. (2.103), and 

£ _1 (n) = £ _1 (n + 1) + 1 (2.106) 

in the case of eq. (2.104). The conditions of eqs. (2.105) and (2.106) are obviously mutually 
exclusive. By repeating the same exercise for L-flows, and by using eq. (2.100), one finally 
arrives at: 

J2 M {n \a',a)= £ <*(£', EJM&^E'.E), (2-107) 
a^ePk s,s'ep; +1 

where A^bxt is given in eq. (2.93), and 
<5(£',£) = J2 uPs(z'~\n)X~\n+l)+a}5(z'~ 1 {n)X~ 1 (n+l)+py (2.108) 

a=-l,l/3=-l,l 

The 5 symbols that appear on the r.h.s. of eq. (2.108) are the ordinary Kronecker delta's. 
Note that only one of the four terms in the sum in eq. (2.108) can be different from zero 
at given £ and £'. These four terms corresponds to those on the r.h.s. of eq. (2.100); their 
signs there are equivalent to the factor a/3 in eq. (2.108). Furthermore, £(£',£) is equal 
to zero when n is not contiguous to n + 1 in either the L- or R-flow. These properties of 
eq. (2.108) ensure that eq. (2.107) holds. Putting all together, one obtains: 

£ Mol^ (£',£), (2-109) 

' n+1 

.Mfe^ (£',£) = r -^— <5(E',E)jp 9fl (z)MW T ( E ' jE ) (2.110) 

i I/O / \ | (knkn+l) K An) /w . [^n&n+l] . >(n) /w | I 



where 



M%> TX - X (Z>, £) = (^ TA (E0|< x ^(S)), (2.111) 
or, in terms of scalar quantities: 

A<S TA ,(S',S) = 4")(E^ r rC7(S',S)^(S J ^). (2.112) 

The consistency between eqs. (2.90) and (2.110) in the sense of eq. (2.29) can be proved 
by direct computation, using eq. (2.43). One can verify that the Kronecker delta's that 
appear on the r.h.s. of the latter equation combine effectively to give <5(£',£) defined in 
eq. (2.108); the quickest way to see this is that of observing that the linear combination of 
the four eikonals in eq. (2.90) can be rewritten as follows: 

£ £ a/3 [£'(£ /_1 (n + 1) + a), £(£^ + 1) +/?)], (2.113) 

a=-l,l/3=-l,l 

i.e. a similar form as eq. (2.108). It is interesting to notice that the collinear limit of an 
eikonal, eq. (2.43), constrains n to be contiguous to n+ 1 in either the R- or the L-flow, but 
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not necessarily in both. However, when n and n + 1 are contiguous only in one of the two 
flows, the linear combination in eq. (2.113) contains two divergent eikonals with opposite 
signs, and therefore is ultimately non-singular. This is the reason why eq. (2.108) forces n 
and n + 1 to be contiguous on both sides of the cut. 

I point out that the treatment of the collinear limit as done above is quite similar to 
that of the soft limit owing to the use of eq. (2.39) rather than of eq. (2.35) as a starting 
point, and in particular to the presence of the Q b operators in the reduced n-gluon matrix 
elements that appear in the former expression. The same technical trick will be adopted 
when treating the case of quark-gluon amplitudes. 

2.3.2 Fixed Born flows 

I now discuss the second of the strategies outlined at the beginning of sect. 2.3. I shall 
make extensive use of the results obtained in sect. 2.3.1, and I shall start from considering 
the soft limit. 

Fixing Born flows is equivalent to the following interpretation of eq. (2.80): at a given 
a, one considers all (n + l)-gluon flows £ = I+(j)cr that arise when performing the sum 
over j. Equation (2.78) then implies that such S's will constitute the set: 

n n 

((a) = |J {l + (a-Hl)) «}^U a} . (2.114) 

i=i j=i 

The discussion that follows eq. (2.80) means that: 

cm n <> 2 ) = if *i^*2, (2.H5) 

U C(*) = K + i ■ (2-116) 

Therefore, the sets C(°") achieve a non-overlapping partition of the space of the (n+l)-gluon 
flows into subsets associated with n-gluon flows. One is thus led to define (n + l)-gluon 
matrix elements at fixed n-gluon flows: 

M {n+1 \a',a) = Yj -M (n+1) (£',£), (2.117) 

sefM £'eC(<x') 

M (n+i) = M (n+1 \a',a), (2.118) 

where eq. (2.118) follows from eqs. (2.115) and (2.116). The derivation carried out in 
sect. 2.3.1 then implies that the soft limit of M^ l+1 \a',a) is equal to M soft'' (o 7 , v) defined 
in eq. (2.69). The colour-linked Born's of eq. (2.70) read, in terms of scalar quantities: 

M<V» = -2A^(a'r[c {l + (a'-\k))a\ I + {a-\l))a) 

-C(l + {o'-\k))o',I-{o-\l))o) 

-C(l-{o'-\k))o',I + {o-\l))o) 

+ C(l-{a'~\k))a', /-((j-^ajjl^V). (2-119) 
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The case of the collinear limit is analogous to that of the soft one. The relevant matrix 
element, Mca^}(a', a), has already been defined, see eq. (2.96). Since the n-gluon flows 
are kept fixed here, the explicit expansion performed in eq. (2.100) is not needed, and one 
can exploit eq. (2.38) to rewrite: 

Mo + llV>) = TT-f— |^(4M (n V>) (2-120) 

+ 2 Q °>- W { [M^"^ ' *> + (Mwi) • " ] ) | ' 

In terms of scalar quantities, A4^ n \a',a) is given in eq. (2.20), and 

M^(a',a)=A^\arC(a',a)Af(a). (2.121) 

On the other hand, the expansion of eq. (2.100) serves to define the analogue of the set 
C(cr), namely: 

Co(<0 = {l + (<r-\n)) a , I-{a~\n)) a] . (2.122) 

At this point, the quantity whose collinear limit is given by eq. (2.120) can be defined in 
the same way as that in eq. (2.117), namely: 

Yl -^ (n+1) (S',S). (2.123) 

However, this is not particularly convenient, since the soft limit of eq. (2.123) is not 
eq. (2.69), which renders it impossible to properly define the subtractions of eq. (1.4). 
On the other hand, one observes that 

Cc(cx) C C(<t), (2.124) 

and that all flows belonging to the (generally non empty) set 

((a) \ C c (a) (2.125) 

correspond to matrix elements which are non-singular in the collinear limit. These two 
facts imply that not only the quantity in eq. (2.123), but also Ai^ n+1 \a' , a) defined in 
eq. (2.117) has the collinear limit given by eq. (2.120), or by its fully equivalent form 
eq. (2.96). One can then conclude that, at fixed a and a', the quantities M. ( - n+1 \cr' , a) , 
{a', a), and Mca^}((T , ,a) have all the properties required to play the same roles 
as the first three matrix elements in eq. (1.4). As far as the soft-collinear counterterm is 
concerned, it is well defined thanks to the consistency of eqs. (2.69) and (2.96) in the sense 
of eq. (2.29). This can be proved exactly as done in sect. 2.1, thanks to eq. (B.14), i.e. to 
colour conservation at fixed Born flows. 
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3. Quark-gluon amplitudes 



I now turn to discussing the case of amplitudes where both quarks and gluons are present. 
The real-emission process will contain q quark-antiquark pairs, and n + 1 gluons: 

— > 2g+(n + l). (3.1) 

The i th parton entering this process will be labelled according to the following conventions 3 : 

—2q < i < —q — 1 — > antiquarks , (3.2) 
—q<i<—l — > quarks, (3.3) 
1 < i < n + 1 — > gluons . (3.4) 

An equal-flavour qq pair will correspond to indices i = —p and i = —p — q; for the moment, 
I shall limit myself to discussing the case of one pair per quark flavour. The quarks 
may be massless of massive; there is no need to distinguish these two cases, since the 
FKS subtraction formulae do not explicitly depend on quark masses (with the obvious 
exception that a massive qq pair cannot induce a collinear singularity - in the following, 
I shall always understand that a quark or an antiquark is massless when is involved in 
a collinear branching). The structure of the underlying Born amplitudes will depend on 
the type of limit considered. When one of the gluons becomes soft, or two gluons become 
collinear, or one quark/antiquark becomes collinear to a gluon, the Born process is as 
follows: 

— ► 2q + n. (3.5) 

In other words, this is the situation in which the FKS parton is a gluon; consistently with 
what was done in sect. 2, such a gluon will have label n + 1. Its FKS sisters will be labelled 
by n, —q, and — 2q for g — > gg, q — > qg, and q — > qg branchings respectively. At the Born 
level, gluon number n + 1 in eq. (3.4) will not appear. The analogues of the labelling of 
eqs. (3.2)-(3.4) will thus be: 

— 2q<i<—q—l — > antiquarks, (3.6) 
—q<i<—l — > quarks, (3.7) 
1 < i < n — ► gluons . (3.8) 

On the other hand, in the case when the FKS parton is a quark, i.e. for a g — > qq branching, 
the underlying Born will be a process: 

— »• 2(g-l) + (n + 2). (3.9) 

The analogues of the labelling of eqs. (3.2)-(3.4) will thus be: 

— 2q + \ <i < —q—\ — > antiquarks , (3.10) 

— q+l<i<—l — > quarks, (3-11) 

l<i<n + 2 — ► gluons. (3.12) 

3 The label i = is not associated with any parton. In the following, I shall not bother to exclude 
explicitly i = when summing over parton labels; this condition will be understood. 
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The FKS parton will have label —q, and its sister label — 2q (which implies that, for this 
case to be non trivial, these quarks must be massless). At the Born level, the gluon with 
label n + 2 will be identified with the one branching into the qq pair. When describing the 
properties of a generic amplitude in the remainder of this section, I shall use the labelling 
of eqs. (3.6)-(3.8) to be definite. 

The colours of gluons will be denoted as in eq. (2.3), whereas in the case of quarks and 
antiquarks: 

{ai}-l_ 2q , <H€{1,...N}. (3.13) 

An amplitude at a given colour configuration, i.e. the analogue of eq. (2.4), will be written 
in the following way: 

A^' n \a. 2q ,...a n )= ]T A({a,},7)I (2 ' ?;n) (7), (3.14) 

where 7 and Tiq;n play the same role as a and P' n in eq. (2.4); namely, they denote a 
flow and the set of all flows relevant to the process of eq. (3.5) respectively. In order to 
determine the precise forms of these quantities, I use again the representation of ref. [1]. In 
particular, one sees that the generic form of the colour structure that multiplies the dual 
amplitude is: 

A ({oi}, 7) = (aM'o+u . . . AM*i)) 

x (\ a «(H+i) . . . A aCT ('2)) 

X . . . 



X 



f\ a "( t q-l+ 1 ) . . . \ a °(tq)\ , (3.15) 



with 

q 

p(j) =max jg- 1, ^S^-p- q,fi(-p- } . (3.16) 

In eq. (3.15), a and fi denote permutations (including cyclic ones) of the first n and q 
integers respectively, a G P n and jjl G P q . The range of the arguments of, and the values 
assumed by, the latter have been extended, so as: 

V(p) = ~V(-P- V) - q, -2q<p<-q-l, (3.17) 
f i(p)€{-q-l,...-2q}. (3.18) 

Finally, the set of q + 1 integers t p 

t = {t ,...t q } , (3.19) 
= t < h < ...t q -i < t q = n, (3.20) 

achieves a partition of the ordered set of the first n integers into q subsets of ordered 
integers, with the p th cell of the partition defined to be: 

(fp_i + 1, t„_i+ 2,... tp-Mp), l<p<q. (3.21) 
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The ensemble of such partitions I shall denote by 

T n \ q . (3.22) 

It should be stressed that the p th cell of the partition, eq. (3.21), may coincide with the 
empty set, which one can formally write by setting 

tp-i = t p . (3.23) 

This case of a zero-length cell corresponds to defining: 

(A^)...AM«) =Sa. p a M (3.24) 

as prescribed in ref. [1]. Note that all cells of t (in the case of amplitudes that feature only 
quarks), or all cells except one (when there is at least one gluon entering the process) can 
have zero length. Equation (3.15) implies the following definition of flow: 

<? 

7 = (J 7p , (3-25) 



P =i 

7 P = 



p ; a(t p -i + 1), . . . a(t p ); /x(-p - qfj . (3.26) 

The physical interpretation of these equations is straightforward. Equations (3.26) repre- 
sents a colour antenna that connects the colour of quark — p with the anticolour of antiquark 
fi{—p—q); attached to this line there are t p — t p -i gluons (therefore, in the case of eq. (3.23), 
there will be no gluons). Then, the integer ^(7), defined in eq. (3.16), is the number of 
times when a quark-antiquark colour line coincides with a flavour line (minus one, when 
there is a maximal coincidence). The flow, eq. (3.25), is the set of all colour lines. By 
convention, the p th colour line begins with quark —p. The construction above implies that 

^2q;n = {Pn,Pq,T n \ q ) . (3.27) 

At this point, one can trivially extend the definitions of colour vectors and the correspond- 
ing amplitudes given in sect. 2. One has just to formally replace 1 . . . n with — 2q . . . n, and 
use the definition of flows given here rather than that of sect. 2. So for example 

\A (2 ^(a^ q , . . . a n )) = A^ n \a_ 2qj ...a n ) |a_ 2 „ . . . a n ) , (3.28) 
|^;n) (7)) = £ A({a i }, 1 )A^ n \ 1 )\a- 2q ,...a n ), (3.29) 

K}?=-2 9 

and so forth. 

The case of amplitudes that feature more than one qq pair per (at least) one flavour 
is essentially identical to what has been discussed so far, and only requires additional 
information on flavour lines, which now are not unique. This information can be given 
e.g. in the following form: 

q 

f= (J (-P,/- P ), (3-30) 
p=i 

{/_!,.../_,} = {-q-l,...-2q}, (3.31) 
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where (— p, f- p ) represents the flavour line that connects quark — p with antiquark 
Clearly, this notation encompasses the equal-flavour case as well, for which /_ p = — p — q 
for all p's. At this point, one can define amplitudes at a fixed flavour configuration /, which 
are gauge invariant, and then sum over all these configurations to obtain the physical 
amplitude. The definition of the colour structure given in eq. (3.15) is unchanged, but 
eq. (3.16) needs be generalized to read: 

<? 

p(j) = max £ 6 n(-p - q)j } . (3.32) 

P =i 

The treatment of equal-flavour amplitudes does not pose any problems, but slightly com- 
plicates the notation. For this reason, in the following I shall deal explicitly only with the 
unequal-flavour case; the equal-flavour case can be easily recovered by adding a sum over 
the flavour configurations given in eq. (3.30). 

3.1 Subtraction of colour-summed matrix elements 

I shall extend here what was done in sect. 2.1 to the case of quark-gluon amplitudes 
squared. The formulae for the relevant soft and collinear limits are of course well known, 
and thus I shall limit myself here to casting them in a form suited to the calculations that 
I shall perform in the rest of this section. As already pointed out before, for quark-gluon 
amplitudes one needs distinguish the cases of the FKS parton being a gluon or a quark. 

I start from the former case. The underlying Born dynamics is therefore that of 
eq. (3.5). The formulae for the soft limit of the amplitude squared are given in eqs. (2.30)- 
(2.32), with only formal changes due to relabelling, and to the fact that self-eikonals need 
not vanish any longer, owing to the possible presence of massive quarks: 



kj 



k,l=-2q 



1 

MgTr +1) = \ 9 2 s £ [Kl]Mtr\ (3-33) 

-(l-<y*l<W). (3-34) 



n+l 



M (2q;n) = _ 2 ^(2,;n)| £ Q\k)Q b (l) \A^) . (3.35) 

b 

The colour operators Q b (k) are defined in eqs. (2.33) and (2.34) in the case when k > 0, 
i.e. when k is a gluon. This definition needs be extended to the case of quarks and anti- 
quarks. We have (see e.g. ref. [11]): 

(a k a n+1 \Q b (k)\c k ) = 6 ban+1 (Q b (k)) , (3.36) 

\ / a k c k 

\ b akCk -q<k<-l & /cisa quark , (3.37) 



(Q\k)) 



akCk 

Q b (k)) = -X b Ckak - 2q < k < -q - 1 k is an antiquark . (3.38) 



With this, one can generalize eq. (2.38) to read: 



Y,Q b (k)Q b (k) = C(k)I, (3.39) 
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C(k) = C A k > , (3.40) 

C(k) = C F £:<0. (3.41) 

The collinear limit can then be cast in the same form as in eq. (2.39): 



Mf q ' n) = (A {2q;n) \ Y J Q\s)Q b {s) \A {2q ' n) ) , (3.43) 

b 

= (A ( ^ n) \ Q b (s)Q b (s) \-A+ q '' n) ) > (3-44) 

b 

where z is defined as in eq. (2.37), with k n — > k s there. As was discussed at the beginning 
of sect. 3, I shall take 

s = n,-q,-2q (3.45) 

as representatives of the cases of the g — >■ gg, q — >■ qg, and q — >■ qg branchings respectively. 
I denoted by X s the identity of the relevant branching parton (i.e., a gluon, a quark, and 
an antiquark) and in eq. (3.42), consistently with eq. (2.42), I have defined: 

Px 3 x s (z) = J^P XsIs (z) , Qx s x t (z) = -^-Qx s x t {z) . (3.46) 

The consistency between eqs. (3.33) and (3.42), in the sense of eq. (2.29), can be proved in 
exactly the same way as was done in sect. 2, thanks to the fact that the colour-conservation 
condition is fulfilled in the case of quark-gluon amplitudes as well (see appendix B): 

n n 

Q b (k)\A^) = => Yl Q b {k)\A (2q >^) = -Q b {s)\A^), (3.47) 

k=-2q k=-2q 

k^s 

and because eq. (2.47) generalizes to read: 

Px 3 xM ^ j-^, QxasW ^ 0, (3.48) 

that is, it also holds in the case of q — > qg and g — S- qg branchings 4 . 

I shall now discuss the case of the g — >■ qq branchings. Here, the underlying Born is the 
process of eq. (3.9), and only the collinear limit is relevant (there are no soft singularities 
when the FKS parton is a quark). The analogues of eqs. (3.42)-(3.44) read: 

(2 9 ;n+l) _ 9 2 s /a , v-rj(2( 9 -l);n+2) 



M ^ k^^ q \ Pq9{Z)Mqq 



+^w*(^^^2!r + 1)! ^)^ (3.49) 



— (2( g -l);n + 2) = ( ^ (2(? _ 1);n+2 ) j ^G^G* ^-l)^)) ? (3 . 50) 

be 

^(%-D;n + 2) = {A( 2 (q -iy,n + 2) l ^Gt c G bc \Af q - iy > n+2) ) , (3.51) 



be 



'For the branchings of quarks and antiquarks, the kernels Q(z) are actually identical to zero [13]. 
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where 

P qg (z) = ^P qg (z) , Q qg *(z) = ^Q qg *(z) . (3.52) 

1 F 1 F 

The (non-hermitian) operators G bc are the analogues of the operators Q b introduced before. 
In colour space, their action is defined as follows: 

(a-2 q a- q \ G bc \a n+2 ) = 8 ba _ q S C a_ 2q Xt_+ 2 a _ 2q ■ (3.53) 

In other words, the operators G bc annihilate gluon number n + 2, and create the qq pair 
with labels — q and — 2q and colours b and c respectively. Note that eq. (3.49) does coincide 
with the usual form of the collinear limit since, as is easy to see from eq. (3.53), one has: 

Y,Gt c G bc = T F I, (3.54) 

be 

which is the analogue of eq. (3.39). 

3.2 Subtraction at fixed colour configurations 

The procedure here is identical to that followed in sect. 2.2, except for trivial changes in 
notation. Its essence is that of inserting into the matrix element squared the projector onto 
a given colour configuration, that in this case reads: 

\a- 2q , ■ ■ .a n+ i)(a- 2q , ■ ■ -a n+ i \ . (3.55) 

After doing that, one will also have to exploit the definitions of the Q h and G bc operators 
given in eqs. (3.36) and (3.53). The quantities whose soft and collinear limits one needs to 
construct are therefore 

M (2*;n+i) {a _ 2qj an+l) = ^;n+i) {a _ 2qj fln+i) | 2 (3 56) 

Its soft limit can be obtained from eq. (3.33) through the procedure described above: 



M { 2q ' n+1) < 



1 n 

'soft (a-2 q , ■ ■ ■ a n+ i) = -g 2 s ^ \k, I A^^' n) (a_ 2 q, . . . a n+1 ) , (3.57) 

k,l=-2q 

M% q;n \a- 2q ,...a n+1 )= (3.58) 
- 2(A&ri\ Q a " +1 (k)\a- 2q , . . . a n+l )(a_ 2q , . . . a n+1 \Q a ^ (I) \A^) . 

which are fully analogous to eqs. (2.51) and (2.52). In terms of scalar quantities, eq. (3.58) 
gets rewritten as follows: 

M ( *r\a- 2q ,...a n+l ) = (3.59) 
- 2 s ^A {2q]n) {a- 2q ...b' k ...at ...a n yQ XkXl (an + i;a k ,ai) b/kbl 



A^ n \a^ 2q ...a k ...h ...an), 



- 25 - 



with 

Ql kIl (an + i;a k , ai ) bc = (Q a ^(k)) bak (Q a ^(l)) aic . (3.60) 

This manifestly coincides with eq. (2.54) when k and I are both gluons. In the case of k 
being a gluon and I being a quark, in SU(3) there are 3-8 2 Q matrices, of dimension 8x3 
(for the case when the quark is on the left of the cut, and the gluon is on the right, the 
matrices have dimensions 3x8). When k and / are both quarks, the are 3 2 -8 Q matrices, 
of dimension 3x3. If quark(s) are replaced by antiquark(s), the number of matrices and 
their dimensions will not change, but their forms will (see eqs. (3.37) and (3.38)). 

The case of the collinear limit is very similar. For g — > gg, q — > qg, and q — > qg 
branchings, from eq. (3.42) one obtains: 

MoT L +1) (a- 2(? , • • • a n+1 ) = ^^\P XsXs {z)M { ^ n \a^ q , . . . a n+1 ) (3.61) 



s ' 



(k s k n+1 )—(2q;n). 



with 



+ Qx s x t (z) K ^±^M™+ '(a- 2q , . . . em+i) 
Mf q ' n \a- 2q , . . . o n+1 ) = (^ (29;n) |Q a ™ +1 (s)|a_ 2? , . . . a n+1 ) 



x <a_ 2g , . . . a n+1 \Q a ^ (s) \A^) , (3.62) 
x <a_ 2g , . . . a n+1 \Q a »+ 1 (s) \A ( + q;n) ) . (3.63) 



M { £$ (a_ 2 „ . . . a n+l ) = (A { l qyn) \ Q a ^(s)\a- 2q , . . . a n+1 



In terms of scalar quantities: 



M ( s 2q;n \a„ 2q ,...a n+1 )= (3.64) 
J2 b' s YQ XsXs (an+i; a s , a s ) b , bs A^ b s ) , 

b' 3 b s 

Mf q Jl\a- 2q ,...a n+1 )= (3.65) 
^^^'"^{aJi^^'J^Qx^^an+^a^a^fe/fe^^^^aJ^,^) . 

Kbs 

The colour matrices that enter this equation are such that 

^2 ^2 Qi s i s (a n +i;a s , a s ) bc = C(s) 6 bc , (3.66) 

dn+l &s 

which extends eq. (2.60), and is the matrix form of eq. (3.39). 

The case of the g — > qq collinear splitting can be derived starting from eq. (3.49). One 
obtains: 

MoT L +1) (a- 2g , • • • a n+1 ) = - 9 \ - \p qg (z)M q f ~ 1);W+2) (a_ 2(? , . . . a n+l ) (3.67) 



2q 



,A {k-qk- 2q ) — (2( g -l);n+2) r . 

+Q qg * [z) M \ — ^—jM qq _ + (a- 2q , . . . a n+1 ) 
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where 

M% q ~ iy ' n+2 \a. 2q , . . . a n+1 ) = (A^~^\ G^^Ja.^, . . . a n+1 ) 

x <a_ 2(? , . . . a n+1 \G a _ qa _ 2q |.4(2(5-i);n+2) ) j (3 68) 

-r-j(2(q-l);n+2), s_ , , (2( 9 -l);n+2) , ^ , . 

x <a_ 2(? , . . . a n+1 \G a _ qa _ 2q \A™ q - 1)rn+2) ) . (3.69) 
In terms of scalar quantities, eqs. (3.68) and (3.69) read: 

M { * {q - iy ' n+2 \a- 2q ,...a n+1 )= (3.70) 
J- A (2(q-i);n+2) ^ 2q+l _ _ ^ a _^ +i _ _ _ a n+1 , b' n+2 yg (a_ q , a- 2q ) b > n+2 b n+2 

x A (2(q ~ iy ' n+2) (a_2 y +i ...X,a_,+i...a n+ i,6 n+2 ), 

(3.71) 

a ™+i,&'„ + 2)*^(a- g ,a-2g)6; +2 fe n+2 

• On+1, ^n+2)j 

with 

G(a- q , a- 2 q)bc = ^ b a- 2q a- q K- q a-2 q ■ (3.72) 

As is clear from eq. (3.72), there are 3 2 Q matrices, of dimensions 8x8 in SU(3). They are 
such that 

^2 ^( a -9> a -2q)bc = T F 5 bc , (3.73) 

which is the matrix form of eq. (3.54). 
3.3 Subtraction at fixed flows 

Although slightly more complicated than in the case of gluon amplitudes owing to the pres- 
ence of quarks, the derivation of the subtraction at fixed flows for quark- gluon amplitudes 
proceeds exactly as that carried out in sect. 2.3. More specifically, I shall use the expres- 
sions of the soft and collinear limits of the colour-summed matrix elements squared, and 
expand them in terms of either the real or the Born flows, using fixed- flows amplitudes and 
the properties of the colour operators Q b and G bc . All the necessary technical ingredients 
are reported in appendix B.2. There, I extend (if necessary) the definitions and derivations 
of appendix B.l, relevant to gluon amplitudes, in a way which renders it straightforward 
the use of arguments by analogy. 

There is however one type of singular limit of the quark-gluon amplitudes that does not 
have an analogue in the case of gluon amplitudes, namely that due to g — > qq branchings. 
I shall deal with these singularities in sect. 3.3.3. 



-^(2(9-1)^+2) 

M qq - + (a-2q, ■ ■ ■ a n+ l) 



b' n+2 b n +2 



(2(g-l);n+2), 



^ -(2(g-l);n+2)/ 

x Ay H (a-2q+i ■ ■ .a^,a- q+ i 
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3.3.1 Fixed real flows 

The quantity whose limits at given real flows I seek to construct is the analogue of that in 
eq. (2.62), namely: 

M {2q ' n+1) (r',r) = (^ 2g;n+1 )(r')|^ {2<?;n+1) (r)) 

= A {2q]n+1 \r'yc(r',r)A i2q ' n+1) (r), (3.74) 

where (see eq. (3.27)): 

r,r' e J" 2?;n+ i = (p n+ i,P ? ,r n+1 |q) , (3.75) 

and the colour-flow matrix elements are defined analogously to eq. (2.21): 

c(r',r)= A {{«}>?)* MM,?). (3.76) 

K}7= + -2 g 

As was done in sect. 2.3.1, one begins by introducing real-emission-level colour vectors 
associated with an underlying Born dynamics: 

|4 2 xT } (r)) = E A (to}- r ) ^ {2q ' ,n) ( r -^) \ a -*> ■ ■ ■ a «+i> • ( 3 - 77 ) 

The Born-level flow IV^- is obtained from the real-level one T by eliminating from it the 
(n + l) th gluon. More precisely, if 



then 



where 



r=(S,M,T)eJ 2 , ;n+1) T = {T ,...T q } , (3.78) 

IWr = (CT, //, t) G J" 2 q;n, (3.79) 

a = (a(l), . . . a{n)) = (S(l), . . .jn^HT, . . . E(n + 1)) , (3.80) 
fi = M , (3.81) 
t = {to,...*,} , (3.82) 

t p = T p Vp s.t. Tp + 1 < S _1 (n + 1), (3.83) 

t p = T p -l otherwise . (3.84) 

Equation (3.80) is the same as in the case of gluon amplitudes - the relative ordering of the 
gluons has nothing to do with their belonging to a given colour line. Eq. (3.80) states the 
obvious but crucial fact that by removing one gluon the structure of the colour lines is not 
affected. Eqs. (3.80)— (3.84) imply that all colour lines in contain exactly the same 

particles and in the same order as those in T. The exception is the colour line derived from 
the colour line f r £ T that contains gluon n+1, since the latter gluon must be removed by 



with 
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definition; the remaining particles on that line, however, have the same relative ordering 
as in T r . 

At this point, as was done in sect. 2.3.1 I start by expressing the soft limit of the matrix 
element squared as a sum over Born-level flows. The analogues of eqs. (2.68)-(2.70) read 
as follows: 

E -^^ +1) (7',7), (3.85) 
\gl E [M]M? n) (7',7), (3-86) 

k,l = -2q 

-2(A {2q ' n) (j')\ E Q\k)Q h (l) \A {2q ' n) (j)) . (3.87) 

b 

The colour-linked Born's of eq. (3.87) can be computed using the results of appendix B.2. 
In particular, from eqs. (B.34)-(B.36) one obtains: 

EQ 6 (oi^ (29;n) (7)) 

6 

EQ 6 (0l-4 (2 " ;ri) (7)> 
b 

Y,Q\i)\A^ n \i)) 
b 

for the cases when I is a gluon (1 < I < n), a quark (— q < I < —1), or an antiquark 
(— 2q < I < — q — 1) respectively. I have denoted by 7 -1 (£) the position of particle / 
in the list that defines the colour flow (see appendix B.2, also for the definitions of the 
operators I±). One can now plug eqs. (3.88)-(3.90) into eq. (3.87) and proceed to the 
explicit expansion of that equation. This calculation can be performed using the very same 
method as in appendix B.2, since the colour structure of the sums over k or I is the same 
as that relevant to colour conservation. More explicitly, one starts by splitting the sum 
over particle labels into sums relevant to single colour lines, and then converts sums over 
particle labels into sums over particle positions. This implies the following manipulations 
(where I consider the sum over I to give a definite example): 

n q ^ q I \ 

E /(o = E E /(o = E /(-p) + E /(Ofei* + m-p - 9)) (3-9i) 

l=~2q p=l le-/ p p=l \ l£y p J 




(3.92) 



where the three terms on the r.h.s.'s of eqs. (3.91) and (3.92) correspond to the quark, 
gluons, and antiquark contributions respectively, relevant to the p th colour line. The ar- 
guments of /() in eq. (3.92) imply that the operators ^±(7 _1 (7(i))) = I±(j) w ih result 
from eqs. (3.88)-(3.90). For each colour line, there will be t p — t p -i + 1 contributions due 



l SOFT 



Ms 2 o^ +1) (y,7) 



^lr } (7',7) 



= |^ n) (/+(7 _1 (0)7)> " l4 2 x 9 T n) (/-(7" 1 (0)7)) , (3-88) 
= |4xT ) (/+(7- 1 (0)7)), (3-89) 
= -|4x 9 T n) U-(7- 1 (/))7)), (3.90) 
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to J_|_ operators (from the gluons and the quark), and t p — t p -i + 1 contributions due to 
i_ operators (from the gluons and the antiquark). Furthermore, the i_ operators can be 
rewritten in terms of i+ operators by using the identity in eq. (B.40). After some trivial 
algebra, fully analogue to that of eqs. (B.39)-(B.42) bar the presence of the eikonal factors 
in eq. (3.86), one obtains: 



. A {2q;n+i) u/ 
- /vt SOFT 



q V t p 

(7,7) = - 9 2 s E E E { [vw.tO')] - [y(o,7(j © i) 

p',p=i i=y _ 1 +iei j'=t p -i+iei 



-/(« © 1), 7(i)] + [V(» © 1), 7(J © 1)] } 



x (4xT ) (/ + «7')l4xT ) (/+(j)7))- (3-93) 



This result is basically identical to that of eq. (2.78), which can be easily understood as 
follows (for the sake of this argument, it is sufficient to consider the case of flows with 
one colour line). In a gluon amplitude, the (n + l) th gluon appears in a given position in 
a real-level flow as the result of the actions of the operators I + and L_ associated with 
the gluon that precedes and follows it in the flow respectively. When the (n + l) th gluon 
occupies the first (last) position, the relevant I + (/_) operator is that associated with the 
last (first) gluon in the flow, as a consequence of the cyclicity of the trace. In the case of 
quark-gluon amplitudes, the colour structure is not cyclic, so the latter argument is not 
valid. However, when the (n + l) th gluon is the left-most (right-most) gluon of the flow, 
there exists a relevant i+ (i_) operator, associated with the quark (antiquark) rather than 
with the last (first) gluon in the flow. It should finally be stressed that the fact that the 
quantities i + 1 and j + 1 which appear in eq. (2.78) are replaced by i © 1 and j ffi 1 in 
eq. (3.93) is only a formal difference. The physical meaning is in fact the same: in both 
cases, these denote the particles to the immediate right of particles i and j. 

Equation (3.93) is in a form suited to transform the sum in eq. (3.85) in a sum over 
real flows. The procedure is identical to that of sect. 2.3.1 and, as in that case, stems from 
the observation that there is a one-to-one correspondence between a real flow T, and a 
pair (7, j), where 7 is a Born flow, and j is an integer whose range is the same as that in 
eq. (3.93). This follows from the existence and uniqueness of the solution of the equation: 

7+(j)7 = r (3.94) 
...(... 7 (j), n + 1, 7 (jffil),...)... 
= ...(... r(j), F(j © 1), F(j © 1 © 1), ...)... , (3.95) 

as was the case for its analogue, eq. (2.82) or eq. (2.83). The use of the © operator in 
eq. (3.95) in place of the ordinary + that appears in eq. (2.83) does not underscore any 
difference between the two cases, but simply allows one to treat with the same notation 
the cases in which j is a gluon or is a quark (note that the range spanned by j in eq. (3.93) 
implies that j is never an antiquark). The existence and uniqueness of the solution of 
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eq. (3.95) can be easily understood e.g. from the relationships that connect real- and Born- 
level flows according to eqs. (3.80)-(3.84). Furthermore, one has to take into account the 
fact that for the two flows j a = (a a , /j, a , t a ) and jb = W>, h) to be different, it is sufficient 
that only one of the conditions a a / o^; \i a / ^ t a / be fulfilled. By solving eq. (3.95) 
one obtains: 



7 = rVtfr, 

j = rV + ijei, 



(3.96) 
(3.97) 



which also imply that: 

7 0') = r(rV + i)ei), 
T(jffii) = r(rV + i)©i). 

In full analogy with eq. (2.86), the arguments above are such that: 

EE E f{-rJ)= E /(r^r-V + i)ei). 

As a consistency check of eq. (3.100), one can again use a counting argument as done in 
the case of gluon amplitudes. Denoting by J\f(n, q) the number of partitions defined in 
eqs. (3.19) and (3.20), one can prove that: 

( n + q- l\ 



(3.98) 
(3.99) 



(3.100) 



N{n,q) 



(3.101) 



The number of terms that appear in the sums on the l.h.s. of eq. (3.100) is 

q\ n\M(n, q) (n + q) , 



(3.102) 



since there are (n + q) different ways of inserting a gluon in a Born flow (n + q is in fact 
the number of terms of the sums over p and j in eq. (3.100)). On the r.h.s. of eq. (3.100) 
the number of terms is instead: 



/! (n + l)\J\f(n + l,q) . 



(3.103) 



Thanks to eq. (3.101), one sees that the numbers in eqs. (3.102) and (3.103) coincide. 

By using the results above, their analogues for the L-flows, and eq. (3.93), one finally 
obtains the soft limit of the matrix element squared at fixed real flows: 



u (2q;n+l) _ 

M 2 o ? PT +1) (r',r) = 



£ Afe +1) (r',r) 



-9 2 s 



{ [r'(r /_1 (n + 1) e 1), r(r- x (n + 1) e i) 



r'(r' (n + i) e 1), r(r _1 (n + 1) © i) 



(3.104) 
(3.105) 



M (2q ' n) 



(r',r) 



r(r (i! + i)0i),r(r> + i)ei) 

nr'" 1 (n + 1) © 1), T(T-\n + 1) © 1)] }m 

(2?;n) (r')|4^ n) (r)). 



+ 



(2q;n) /-p/ 



(r',r), 

(3.106) 
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The form of eq. (3.105) is fully analogous to that of eq. (2.90). In terms of scalar quantities: 

m^\t', r) = A^Xr'^y c(V, r) A^ir^) . (3.107) 

The derivation of the soft limit of the matrix elements given above shows the complete 
similarity between the cases of gluon and quark-gluon amplitudes. The collinear limits 
relevant to the branchings g — > gg, q — >■ qg, and q — > qg are also fully analogous, as one can 
easily understand by comparing eqs. (2.39)-(2.41) with eqs. (3.42)-(3.44). In other words, 
in the case of quark-gluon amplitudes one can follow the procedure that in sect. 2.3.1 has 
led me from eq. (2.95) to eqs. (2.109)-(2.111). I refrain from doing that explicitly, and 
limit myself to present the final result: 

MoT L +1) = £ MoT L + V,r), 

r,r'6J"2(j;n+l 

M {2q '' n+l) <r' r\ - 



l COLL 



K s' K n+l I 



+\Q^M ( %^mSUF,t) + ^^M^Ur ,r) 




where 



MS(r',r) = (A^(T')\A^S(T)), (3.110) 

and 

J(r',r)= Yl a/35(r ,_1 (s),r /_1 (n + l)©a) 5(T nl (s),T nl {n + l)(B^y (3.111) 

a=-l,10=-l,l 

In terms of scalar quantities: 

^S(r',r) = If ;n) (r^-r c(r',r)lf^(r^) . (3.112) 

The very close correspondence between eqs. (3.105) and (3.109), and their counterparts 
in the case of gluon amplitudes, implies that the consistency (in the sense of eq. (2.29)) 
between the soft and collinear limits of quark-gluon matrix elements can be proved in 
exactly the same way as explained at the end of sect. 2.3.1. 

One comment is in order as far as the form of eq. (3.111) is concerned. In the case of 
a 9 ~^ 99 branching, i.e. when s = n, the situation is strictly identical to that explicitly 
presented in eqs. (2.100)-(2.102). However, when one considers a q — > qg branching, i.e. 
when s = —q, the analogue of eq. (2.100) contains only one term, because of eq. (3.89). 
Therefore, the relevant real flows, which correspond to those in eqs. (2.101) and (2.102), 
are: 

I + irT\-q))i, I + {T\-q))l — > (-q;n+l,...). (3.113) 

These correspond to the contribution obtained when a = —1 and (3 = —1 in eq. (3.111). 
However, it is easy to see that if either a = 1 or f3 = 1 in eq. (3.111), then <5(r',r) = 
if s = —q. In fact, if 6(T', T) = 1 these cases would correspond to real flows where the 



- 32 - 



(n + l) th gluon occupied the position at the immediate left of quark —q, which is impossible. 
Similar arguments apply to the q — > qg branching, i.e. when s = —2q. The bottom line is 
that eq. (3.111) can indeed be used to describe all branchings considered here, as claimed 
before. 

3.3.2 Fixed Born flows 

The results relevant to fixing the Born flows can be straightforwardly obtained from 
sect. 3.3.1. Once again, the analogy with the case of gluon amplitudes is very close. 

The soft limit of the matrix element has already been presented in eqs. (3.85)-(3.87). 
The real-emission-level quantity whose soft limit is M^o^ +1 \'-f' , ^/) w iU be a linear combi- 
nation of the fixed-real-flow matrix elements Ai^ 2q ' n+1 ^ (V ,T) defined in eq. (3.74): 

.M (2<z;n+1) (y,7) = Y, Yl -M (29;n+1) (r',r), (3.H4) 
re£(7) r'e£( 7 ') 

for suitable sets of real flows £(7) and £(7'), to be defined later. In terms of scalar quantities, 
the colour-linked Born's of eq. (3.87) read: 





-<w 




- <w 




- <w 


Sqx k )(l 


- <w 






(3.115) 



xA i2q > n) {~f). 

The (1 — 5) terms that multiply the colour-flow matrix elements in the r.h.s. of eq. (3.115) 
exclude the contributions of the 7_ and 7 + operators in the case when the corresponding 
particle is a quark or an antiquark respectively, as dictated by eqs. (3.89) and (3.90). 

Equation (3.114) is the analogue of eq. (2.117), and the set £(7) is the analogue of 
its counterpart £(cr) in the case of gluon amplitudes, defined in eq. (2.114). Since ((a) is 
the set of all real flows that can be obtained by acting on a given Born flow a with the 
operators 7± relevant to the computation of the colour-linked Born's, one can obtain £(7) 
precisely in the same way. Hence, from eqs. (3.87)-(3.90), one gets: 

n — 1 — q— 1 

e(7) = U{ / +^~ 1 «)7,/-(7" 1 (0)7} U {^+(7- 1 (0)7} [J {^-(7- 1 (0)7} 

1=1 l=-q l=-2q 

= U U { J +0')7}, (3.H6) 
p=i j=t p _i+iei 

where the last form follows from eq. (B.40), but could equally well be deduced directly 
from eqs (3.88)-(3.93). The derivation presented in sect. 3.3.1 implies that: 

£(7i)n^ 2 ) = if 71^72, (3.117) 

(J £(7) = ^;n+l , (3.118) 
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as in eqs. (2.115) and (2.116). Hence: 

M (2q;n+i) = £ M (2q ' n+1 \j',j), (3.119) 

which confirms that indeed the quantity defined eq. (3.114) has the property that one 
expects. 

In the collinear limit, the matrix elements at fixed Born flows can be obtained directly 
from eqs. (3.42)-(3.44), by using the identity of eq. (3.39): 

MoT L +1) (7',7) = TT-f— K* s (^ (29;n) (7',7) (3.120) 

K s' K n+l I 

l r , / J (KK+l) k A2q;n), / x . [^sK+l] K A2q;n), I x\ \ 

+ 2 Q ^ Z \jk^] M -- (7 ' 7)+ (iUW)^- (7 ' 7) J)' 
with the reduced matrix elements that have the usual definitions: 

M^ n \ 7 ',j) = (^ (29;n) ( 7 ')l^ (29;n) (7)) 

= ^ 2<? ' n )( 7 , ) ,t C(7 / ,7)^ {2 ' ?;n) (7) , (3.121) 
M 2 fV,7) = {A^ n \j')\A^ n \ 7 )) 

= A ( ^ n \ 7 'yW, 7 ) AW ( 7 ). (3.122) 

The very same arguments as in the case of gluon amplitudes can now be applied. Namely, 
one can construct the sets 

Ul) = {/ + (7~V))7, /-(7~V))7} , {/+(7- 1 (-^))} , {i-^-H-m} , 

(3.123) 

relevant to the g — > gg, q — > qg, and q — >■ qg branchings respectively. However, one finds 
again that 

Cc(7) C £( 7 ), (3-124) 
for the three cases, and that all flows belonging to the sets 

£(7) \ Cc(7) (3-125) 

do not induce collinear singularities. Therefore, the matrix elements defined in eq. (3.114) 
have the collinear limits given by eq. (3.120), and can thus be used with the latter and 
with those of eq. (3.86) for the subtractions of eq. (1.4). 

3.3.3 The g — ^ qq branching at fixed flows 

As was the case for the other collinear branchings discussed so far, the starting point is 
the expression of the colour-summed collinear limits (here, eqs. (3.49)— (3.51)), where one 
writes the scattering amplitudes using their representations in terms of Born colour flows. 
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Therefore: 



- /V/l COLL 



£ M 2 oT L +1) (7',7), 



7,7'eJ r 2 ( q -i);„ + 2 



A^V,7) = y^f— ^ ^(^)^ 1);n+2) (7',7) 



k— q 'k—2q 



(3.126) 
(3.127) 



+ 7,Qqg*( z ) 



(k- q k- 2 q)_j^(2(q-l);n+2)^ l ^ , [k- q k- 2 q] -yjpfa-l);^) ^/ 



2^ s I [fc_,fc_ 2 ,]"'«' 



'(Y,7) + 



'(Y.7) 



where 



^ (9 " 1);n+2) (V,7) = (^-^(Y)! £ G £Ac |^P(«-i);n + 2 )(7)) j (3 _ 128) 



6c 



(2(g-l);n+2)^, 



(Y,7) = <4 % ~ 1);n+ V)l E G tcG bc l4 2(9 - 1); " +2) ( 7 )). (3.129) 



(2(g-l);n+2), 



be 



Hence (using the same arguments as in eq. (2.74)), the relevant quantities to compute are 
the colour vectors: 

G hc \A {2{q - 1) ' n+2) { 1 )) (3.130) 

(or their bra counterparts), which are the analogues of those in eqs. (3.88)-(3.90) or 
eq. (2.73). By using the definition of the operator Gb c given in eq. (3.53), one sees that the 
colour part of eq. (3.130) involves the computation of: 



— ' V / a_ r a 

O-n+2 



_ ( + . . . / \ a <r(a-l(n+2)-l) j M a CT ( CT -l(n+2) + l) _ _ _ \°-a(t r ) 

2 V / a_ r c V 

— he (A a<T(t '- 1+1) . . . A aCT ( CT_1 ("+ 2 )- 1 )A aCT ( CT_1 ("+ 2 )+ 1 ) ...X^rA 



V(-i — ?) 



-«) 



(3.131) 



where, as the notation suggests, r is the colour line to which gluon n + 2 belongs, and the 
r.h.s. of eq. (3.131) has been computed using eq. (A. 5). Therefore, for any Born flow 



7 £ - 77 2(g-l);n+2 



9-1 

7 = U > ' 
P =i 



(3.132) 



eq. (3.131) suggests to define two operators as follows: 



n + 2 G 7 r , J7 



(J (tf 7r ) , (3.133) 
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where by construction J7, if 7 G i^n+i, and 

J 7r = (-r;a(t r _ 1 + l),...a(^ 1 (n + 2)-l);-2( 7 ) J 

( -g ; a(a-\n + 2) + 1), . . . <r(t r ); /i(-r - g)) , (3.134) 

i^7r = (-r;a(t r _i + l),...a(a" 1 (n + 2)-l), 

C r( ( j- 1 (n + 2) + l),...a(t r );M-r--g)) IJ (-g;-2g). (3.135) 

In other words, the operator J splits in two the colour line to which gluon n + 2 belongs. 
The gluons to the left of gluon n + 2 now belong to a colour line in which the colour of 
the quark of the original line forms an antenna with the anticolour of antiquark —2q (that 
emerges from the g — >■ qq branching) , whereas the gluons to the right of gluon n + 2 now 
belong to a colour line in which the colour of quark — q (that emerges from the g — > qq 
branching) forms an antenna with the anticolour of the antiquark of the original line. On 
the other hand, operator if simply removes gluon n + 2 from the original colour line, and 
creates the colour line that connects — q with — 2q; the latter line does not contain any 
gluons. With the definitions above, one arrives at: 

]TG bc |^-i);^) (7)) = i|4 2 ir 1);n+2) (J7,7)> - ^|4 2 ir 1);n+2) (^7,7)) , (3-136) 

be 

where the definition of eq. (3.77) has been extended: 

|4x ( r 1);n+2) (r,7)>= E A({a0,r)^-^ +2 )(7)|a_ 2(? ,...a n+1 ). (3.137) 

for any T G F 2q - n +i , 7 G ^2( g -i);n+2 • (3.138) 

Note that the factor l/N that appears in the second term on the r.h.s. of eq. (3.131) is 
indeed contained in the second term on the r.h.s. of eq. (3.136) owing to the definition of 
K, eq. (3.135). In fact, the latter equation shows that if 7 features a colour line (—q; — 2q) 
which is also a flavour line, and therefore: 

P (K 1 ) = p( 1 ) + l. (3.139) 

At this points, one plugs eq. (3.136) into eqs. (3.128) and (3.129). Then, the sums over 
Born flows in eq. (3.126) can be turned into sums over real flows using the same technique 
as was employed in eq. (2.80). To be definite, I consider the case of R-flows, that of L-flows 
being identical. One has 

E (/(J7, 7) +5(^7, 7)) = 

7e^2(( 7 -l);n + 2 

E E (5 r ' J 7)iV7, 7) + S(T, Kj)g(K 7 , 7)) , (3.140) 

76^2(9- l);n+2 reJ2q;n+l 
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where I have used the identities: 



1= E 5 ( r > J 7), (3.141) 

reJ-2ij:n + l 

1= E <H r ^7), (3.142) 



which hold because, for a given 7, the flows J7 and exist and are unique. The 5 
functions can then be used to get rid of the sum over 7, which requires fixing T, and 
solving: 

r = J 7 , (3.143) 
T = K~f, (3.144) 

for 7 (these equations are the analogues of eq. (2.82) or eq. (3.94)). I start by considering 
eq. (3.143). As one can deduce from eq. (3.134), eq. (3.143) has a solution if and only if T 
belongs to a subset T^n+i °f se * °f rea ^ fl° ws > defined as follows: 

r=(<7, M ,t)€ 4 J q l +1 CT 2q -n + i /i(-2?)^-2g. (3.145) 

From eq. (3.134), it is also obvious that when such a solution exists, it is also unique, and 
I shall denote it by 7 = J _1 r. Therefore: 

E /( J 7,7)^ E E *(T^7)/(^7,7)= E /( r '^ lr )- 

(3.146) 

Let me now turn to eq. (3.144) which, as can be seen from eq. (3.135), has a solution if 
and only if T belongs to a subset F^q-n+i °f se * °f rea ^ fl ows > defined as follows: 

T = {e,n,t) € C JF 29;n+1 <=► M (-2g) = -2(7 and = t q = n+l. (3.147) 

At variance with the case of the operator J, however, the solution of eq. (3.144) is not 
unique. In fact, the operator K just removes gluon n + 2 from the first q—1 colour lines, 
replacing it with the colour line (—q; — 2q). Hence, the Born flows in which all quarks and 
gluons have the same relative positions, up to that of gluon n + 2, will result in the same 
real flow if 7. It is easy to convince oneself that this happens exactly (n + q) times. I shall 
therefore define the multi-valued inverse of the operator K, which will give the solutions to 
eq. (3.144), as follows: 

^'re^-i)^ for any 1 < i < n + q , r 6 j£ +1 . (3.148) 

The explicit form of K~ l T can be worked out from eq. (3.135); it corresponds to inserting 
gluon n + 2 in all possible ways in T, except in the colour line (— q, —2q), which is removed 5 . 



3 This implies that K 1 could be written in terms of the operators 7+ defined for gluon n + 2, if need be. 
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Putting all this together, one obtains: 

n+q 

g(Kj,j)= E S(t,k 1 ) 9 (k 1 , 1 ) = E E^^r)- 

(3.149) 

Counting arguments can again be given as a consistency check of the results obtained here. 
The number of elements in the sets of flows relevant to the derivation above are: 

# (^2(,-i);n+2) = (? - 1)! (n + 2)UV(n + 2, g - 1) , (3.150) 
#(^+1) =(g-l)(g-l)!(n + l)!JV(n + l, 9 ), (3.151) 
= (g-l)!(n + l)U\f(n + (3.152) 

with N given in eq. (3.101). By direct computation, one shows that 

# (^2( ? -l);n + 2) = # (4'ql+l) , (3-153) 

# (^2( ? -l);n+2) = # (^S+l) (" + ?)■ (3-154) 

Hence, the same number of terms appear on the two sides of eqs. (3.146) and (3.149). 
Equations (3.126)-(3.129) can now be expressed in terms of real flows: 

Afe +1) = E MoT L +1) (r',r), (3.155) 

r r'aT^ I I 

M^Ll +1 \T>,r) = ^ & k 2 ^ P qg (z)M^ n+2) (T',T) (3.156) 

_j_ 1 A ( 7 \ ( ( k -q k -2q) . ,(2( g -l);n+2) rr / , [fc-qfc-2g] ■ . (2(g-l);n+2) /w px A \ 
+ 2^fl*W ^ [fe _^ ^ -MEXT.AB.-+ (1 (L?L2 ^ E XI,AB, + - 7 J ' 

where, in eq. (3.156): 

A,B G {J, K} , r' € J" 2 ( ^ +1 , T G . (3.157) 

The reduced matrix elements squared are defined as follows: 

Mxtn );n+2) (r',r) = (^'^^(r'^^rOl^ir 1 ^^^,^^)), (3.158) 

n+q 

Mxt:K );n+2) (r',r) = -^(4x { r 1);n+2) (^^J- 1 ^0|4 2 x ( ^ 1);n+2) (^,^7 1 ^)), (3.159) 

i=i 

n+q 

-Mgfe:™ );B+2) (r',r) = -X;<^T- 1):n+2) (r / ,^r 1 r , )l^r 1):n+2) (r, J^r)), (s.ieo) 

i=i 

n+q n+<? 

M&S^V.r) = ££<4 2 x ( r 1);n+2) (r^^^^ (3.161) 

1=1 j=l 
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and similarly for the azimuthal terms -M^t ab'ax~ 2 ^> wm cli as usual only require the use of 
amplitudes at given polarizations. The matrix elements of eqs. (3.158)-(3.161) all feature 
the same colour structure, as is implied by eq. (3.137). This becomes evident if one expresses 
them in terms of scalar quantities: 

M 2 x ( tiJ );n+2) ( r '> r ) = AMi-^+V (J^T'y C{V, T) l(2(?-i);n+2) ( j-i T) j (3 162) 

(n+q 
^^(2(,-l);n+2) (K -l r) 

(3.163) 



(3.164) 



M (2(q-l);n+2) /r , 



1 EXT.KK 



(r',r) = y^a^-V'^Xk^t')* c(t',t) 



'n+q \ 
^^(2( g -l);n+2) (jK --l r) \ _ (3 J65) 

As is implied by eq. (3.155), a closed flow (r',T) will not induce a singularity if either its 
L-flow or its R-flow (or both) does not belong to the set j|^ n+1 U-^io-n+i- 

Collinear singularities due to g — > qq branchings can obviously occur when the under- 
lying Born amplitude is a pure-gluon one, which corresponds to q = 1. The notation used 
so far is not suited to describe such a case, the Born amplitudes having been written as 
quark-gluon ones. It is however not difficult to extend what was done before to q = 1, since 
one can formally understand 

|^(2( 9 -l);n+2) (7)) ^ ^(0^+2)^ = 1^+2)^ ^gg) 

with a being the only non-trivial part of the flow 7 in this case, namely a permutation 
of gluon labels (one must also understand that the permutations are restricted to be non- 
cyclic only, when q = 1 is considered). Equation (3.130) is then still relevant, but its colour 
part is not given by eq. (3.131), but rather by: 

^ \^ +2 Tr(^X a "W . . . AV^ 1 ^ 2 )- 1 ) X an + 2 \ a <r(<r- 1 ("+V+l) . . . A<M"+2)^ = 
- ( \ a <r(<T- 1 (n+2) + l) \ a a(n+2) \ a a(l) A"^" 1 ("+ 2 )" 1 ) ^ 

2 V "' ••• j bc 

_ -L^Tr^A"^ 1 ) ...A a -( CT_1 ("+ 2 )- 1 )A aCT ( CT_1 ("+ 2 )+ 1 ) ...AM»+2)^ . (3.167) 

The second term on the r.h.s. of eq. (3.167) does not belong to the set of real flows; hence, 
it must give a contribution to the final result equal to zero. That this is indeed the case can 
be proved by direct computation. One starts by extending the definitions of the operators 
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J and K to the case q = 1: 

J7 = (-g;a(^ 1 (n + 2) + l),... CT (n + 2), 

(7(1), . . . aia'^n + 2) - 1); -2g) , (3.168) 
K 7 = (a(l),...(7(<7- 1 (n + 2)-l) ) 

<7(<T _1 (n + 2) + 1), . . . a(n + 2)) [j ( -g ; -2g) . (3.169) 

With these, one proceeds exactly as done before. The only difference is that now 

= ^n+l ■ (3.170) 

(K) 

As far as F^q-n+i ls concerned, one simply defines it as the set of flows that have the same 
form as the one in eq. (3.169). By doing that, ones arrives again at eqs. (3.155)-(3.161). 
However, the following reduced matrix elements are identically equal to zero: 

A<iS(r',r) = MS(r',r) = MiS(r',r) = o, (3.171) 

since they contain one or both of the linear combinations 

n+1 n+1 

Y^A^^Xk^Y') = Y J ^ n+2) (Kj 1 T) = 0, (3.172) 

i=i j=i 

which are equal to zero owing to the dual Ward identity [1]. This proves the fact that the 
second term on the r.h.s. of eq. (3.167) indeed does not contribute to the final result. The 
collinear limit for the g -+ qq branching in the case of a pure-gluon Born is therefore still 
given by eqs. (3.155) and (3.156), with the condition that T^-n+i ^ e ec i ua l to the empty 
set, and taking eq. (3.171) into account. 

In analogy with what was done previously, I should now discuss the case of the g — >■ qq 
branching at fixed Born flows. However, it does not appear to be possible to formulate it 
in a gauge-invariant manner. This is due to the fact that the operator K does not have a 
single-valued inverse. If one must keep the Born flows K~ l T separate, rather than summing 
them as done in eqs. (3.159)-(3.161), one must also "split" the contributions to the real 
flow r into (n + q) components, in order for each of them to have a limit proportional 
to ^K^T). This may be possible by considering individual Feynman diagrams, at the 
price of violating gauge invariance. It does not seem justified to do so since, contrary to 
the cases of g — > gg and q — >■ qg branchings, the Altarelli-Parisi kernel associated with the 
g — > qq branching is O(N ), which implies that for this branching the two formulations at 
fixed Born or real flows are strictly equivalent from the colour point of view. The above 
discussion does not apply to the case q = 1, i.e. when one has a pure-gluon Born since, as 
was shown before, when q = 1 the K operator does not contribute to the result. But it is 
obvious that in such a case fixing Born flows is identical to fixing real flows. 
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4. Summary of subtraction formulae 

Since the detailed derivations carried out in sects. 2 and 3 can obscure the final results, in 
this section I collect the references to the formulae needed for the implementation of FKS 
subtraction at fixed colour configurations or flows. For the colour-summed real-emission 
matrix element squared M. and its limit A4 L (be it soft, collinear, or soft-collinear) : 

M^M h , L = SOFT, COLL, SC, (4.1) 

I have used the decompositions: 

M = Y,M{g), M L = Y,M L (g), (4.2) 

9 9 

with 

M(g)^M L (g) V 5 . (4.3) 

I have considered three cases, which I list here using the notations relevant to gluon and 
to quark-gluon amplitudes respectively. As discussed in the text, the soft-collinear limit 
A4 sc (g) can be trivially obtained by computing either the collinear limit of M SOFT (g), or 
the soft limit of -McoliXsO; hence, it will not be considered in what follows. 

• Fixed colour configurations: 

g = (a 1 ,...a n+1 ) , g = (a_ 2(? , . . . a n+1 ) . (4.4) 

The matrix elements M(g) are defined in eq. (2.48) for gluon amplitudes, and in 
eq. (3.56) for quark-gluon amplitudes. 

Soft limits Ai SOFT (g). For gluon amplitudes: eqs. (2.51) and (2.53). For quark-gluon 
amplitudes: eqs. (3.57) and (3.59). 

Collinear limits M C oll{9)- For gluon amplitudes, eqs. (2.55), (2.58), and (2.59). For 
quark-gluon amplitudes, and g — > gg, q — > qg, or q — > qg branchings: eqs. (3.61), 
(3.64), and (3.65). For quark-gluon amplitudes, and g — > qq branchings: eqs. (3.67), 
(3.70), and (3.71). 

• Fixed real flows: 

5 = (£',£), g={T',T). (4.5) 

The matrix elements M{g) are defined in eq. (2.62) for gluon amplitudes, and in 
eq. (3.74) for quark-gluon amplitudes. The definitions of flows in the two cases are 
given at the beginning of sect. 2 and of sect. 3 respectively. 

Soft limits .Msoft^)- For gluon amplitudes: eqs. (2.90) and (2.93). For quark-gluon 
amplitudes: eqs. (3.105) and (3.107). 

Collinear limits M. C o^{g)- For gluon amplitudes, eqs. (2.110), (2.93), and (2.112). 
For quark-gluon amplitudes, and g — > gg, q — > qg, or q — >■ qg branchings: eqs. (3.109), 
(3.107), and (3.112). For quark-gluon amplitudes, and g — >■ qq branchings: 
eqs. (3.156) and (3.162)-(3.165). 
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• Fixed Born flows: 

g = (a', a) , g = (7', 7) 



(4.6) 



The matrix elements M{g) relevant to this case are linear combinations of those 
defined at fixed real flows. They are given in eqs. (2.117) and (2.114) for gluon 
amplitudes, and in eqs. (3.114) and (3.116) for quark-gluon amplitudes. 
Soft limits M SOFT (g) . For gluon amplitudes: eqs. (2.69) and (2.119). For quark-gluon 
amplitudes: eqs. (3.86) and (3.115). 

Collinear limits M. C o^{g)- For gluon amplitudes, eqs. (2.120), (2.20), and (2.121). 
For quark-gluon amplitudes, and g — > gg, q —> qg, or q — > qg branchings: eqs. (3.120)- 
(3.122). For quark-gluon amplitudes, and g — > qq branchings: this case cannot be 
treated at fixed Born flows. 

5. Born-like contributions 

I have so far discussed the case of the real-emission matrix elements and of their local 
subtraction terms. NLO cross sections also receive other contributions, which have a 
Born-like kinematics. These contributions are due to the Born proper, to the one-loop 
corrections, and to the analytically-integrated subtraction terms. In this section, I shall 
discuss the treatment of the latter when one fixes the colour configurations or the colour 
flows. 

In FKS one defines two types of Born-like contributions, which arise from the analytical 
integration of either the soft or the collinear counterterms. The one of soft origin reads as 
follows, up to overall trivial factors (see e.g. ref. [11]): 



where I have used the notation of the quark-gluon amplitudes case for generality. Equa- 
tion (5.1) coincides with eq. (3.33), except for the fact that the eikonal factors [k,l] in 
the latter equation have been replaced in eq. (5.1) by the finite parts of their integrals, 
denoted by S^i- The eikonals do not play any role in the manipulations carried out in 
sects. 2 and 3. Hence, one can just use the results for the soft limits of the matrix elements 
presented before and summarized in sect. 4, and simply replace [k,l] with £m there. 

The Born-like contribution of collinear origin is proportional to the Born matrix ele- 
ment squared: 



n 




(2q;n) 
kl 



(5.1) 



k,l=-2q 



(5.2) 



where 




k=-2q . 



n 



7 (X fc )-2C(J fc )log 
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and 



i{g) = 


— c A - 
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-\t,n,, 


(5.4) 
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3 c 
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i{g) = 


(67 
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The definitions of the various quantities that appear in eq. (5.3) are irrelevant here; the 
interested reader may find them in ref. [11]. What matters is the definition of the colour 
factors, given in eqs. (5.4)— (5.7). These basically arise from the integration of the Altarelli- 
Parisi kernels: P qq for 7(g) and j'(q) (q — S- qg and q — ^ qg branchings); P gg plus P qg for 7(g) 
and 7' (g) (g — > gg and g — >■ qq branchings). The sum of P gg and P qg is necessary because 
7(5) and 'j'(g) are inclusive properties of gluons, and therefore all possible branchings must 
be taken into account. However, these can still be told apart in eqs. (5.4) and (5.6), thanks 
to the different colour factors {C A versus T F ). This implies that, in eq. (5.3), one can 
unambiguously associate the various terms with either a q — > qg or a q — > qg branching 
(C F terms), or a g — > gg branching (C A terms), or a g — > qq branching (T F terms). For 
each of these branching types, one can then repeat what was done in sects. 2 and 3 in the 
case of collinear limits. This implies starting from re-writing eq. (5.2) as follows: 

— 1 n n 

Q M (2q;n) = ^ C F . . M (2q ' n) + C a( ■ ■ ■) M (2q ' n) + Y, T F ( ■ ■ ■) M {2q ' n) , (5.8) 
k=-2q ' k=l k=l 

where the terms in brackets can be easily worked out from eqs. (5.3)— (5.7). Their forms 
are irrelevant here, except for the fact that they do not contain any colour factors. At this 
point, one exploits eq. (3.39) in the first two terms of eq. (5.8), to replace M^ 2q ' ,n ' > with 
Mf q ' n) defined in eq. (3.43) (with s — > k there). Analogously, one exploits eq. (3.54) in the 
last term of eq. (5.8), to replace M^ q ^ with M q 2 q (q ~ iy ' n+2) defined in eq. (3.50). After these 
replacements, the matrix elements in eq. (5.8) are in the same form as those that appear 
in the colour-summed expressions of the collinear limits. Hence, the same manipulations 
as in sects. 2 and 3 can be carried out here; the final results can be directly obtained from 
the expressions of the collinear limits at fixed colour configurations or real flows, by simply 
replacing the Altarelli-Parisi kernels and prefactors there with the expressions in round 
brackets that appear in eq. (5.8). Note that this procedure must not be carried out if one 
is interested in fixing Born flows: the original expression, eq. (5.2), is already suited to 
that, the only change being the formal replacement: 

M (2q ^ — > M {2q ' n) (V,7). (5.9) 

6. Discussion 

Processes with large particle multiplicities have a colour algebra so involved that its direct 
computation is impossible, and one must use alternative methods, such as Monte Carlo 



(MC) ones. While it is common to sum over colour configurations with MC techniques, 
there is no reason of principle that prevents one from doing an MC sum over colour flows - 
such a strategy is indeed being considered in the new version of MadGraph/MadEvent [31]. 
At the tree level, matrix elements at fixed colour configurations (e.g. eq. (2.18)) have the 
advantage over those at fixed flows (e.g. eq. (2.20)) of being positive definite, while the 
latter are in general complex numbers. However, this advantage is lost beyond the leading 
order in perturbation theory, since subtractions such as those of eq. (1.4) will always be 
involved, thus implying results of either sign. Also, in an actual computation one will 
consider the sum M(a',a) + M(a,a'), which is a real number, since physical observables 
must not depend on the distinction between L- and R-flows. Furthermore, at any order 
in perturbation theory, including tree level, matrix elements at fixed colour flows have an 
immediate interpretation in terms of large- N expansion, which is rather indirect at fixed 
colour configurations. In particular, it is not necessary to explicitly compute the colour- 
flow matrix of eq. (2.21) (or eq. (3.76)) to determine the largest possible power of N in 
any of its elements - it is sufficient to count the number of colour loops determined by the 
closed flow (a', a) (or (7', 7)), which is basically instantaneous if performed by a computer. 
This fact not only paves the way to a systematic organization of the computation in terms 
of increasing powers of 1/N, but also suggests a way to save computing time, since terms 
with large 1/N powers may be computed with a relatively low statistics. 

The results presented in this paper allow one to follow either of the strategies discussed 
above in the context of the computation of NLO observables with FKS subtraction. The 
approach that uses fixed colour configurations bases its efficiency on the fast calculation 
of colour-dressed amplitudes, at both the real-emission and Born levels. It should be 
stressed that, for a given real matrix element at fixed colours, such e.g. that of eq. (2.48) 
for gluon amplitudes, one must compute several Born-level colour-dressed amplitudes - 
see eq. (2.53) for the soft limit, and eqs. (2.58) and (2.59) for the collinear limit. It is 
clear that, for the sake of numerical stability, the sums over b' k and bi in eq. (2.53), and 
those over b' n and b n in eqs. (2.58) and (2.59), must be performed exactly, and not with 
MC methods. Furthermore, in the case of the soft limits the sums over k and / that 
appear in eq. (2.51) must also be taken into account. The bottom line is that, if one 
performs these sums blindly, a non-negligible complexity creeps back into the game. For a 
fast computation of the counterterms, then, an algorithm is essential that pre-determines 
which colour configurations are associated with tree-level amplitudes that are equal to zero. 
It is probably also important to compute the colour matrices Qx k Xi and Q once and for all, 
since many of their elements will be equal to zero, an information which can be effectively 
used to improve the efficiency of the algorithm. These issues are beyond the scope of this 
paper, and will not be discussed any further. 

Turning to the fixed-flow approach, the general arguments related to the 1/N expansion 
apply to both the fixed-real and the fixed-Born flow cases. I also stress that, in order to 
derive the formulae for the fixed-flow subtraction terms, I have started from the known 
limiting behaviours of the squared matrix elements, the same that I have used in the case 
of fixed colour configurations. In other words, I did not employ the limits of colour-ordered 
amplitudes. While the latter may have provided a more straightforward derivation for fixed 
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real flows, the approach followed here allows the treatment of the three schemes considered 
(fixed colour configurations, Born flows, and real flows) in a common language. The fact 
that the results I have arrived at for the matrix element limits at fixed real flows can also 
be derived using the limits of dual amplitudes as a starting point constitutes a partial cross 
check of the procedure adopted in this paper. 

When fixing the Born flows, the building blocks of the counterterms are the same ones 
as those used in the automatic implementation of the FKS subtraction achieved by Mad- 
FKS [11]. Furthermore, in this scheme the simultaneous integration of the real-emission 
and of the Born-like and Born contributions can be performed in exactly the same way 
as in ref. [11], where it is discussed in detail. On the other hand, from the point of 
view of computational complexity, for any pair of Born-level dual amplitudes that deter- 
mine the counterterms, one needs to evaluate n 2 (in the case of gluon amplitudes - see 
eqs. (2.117) and (2.114)), or (n+q) 2 (in the case of quark-gluon amplitudes - see eqs. (3.114) 
and (3.116)), real-emission dual amplitudes. This has an obvious physical interpretation: 
the sets of real flows ((a) and £(7) effectively achieve a block decomposition of the colour- 
flow matrices C(S',S) and C(T',T) respectively. The choice of a Born-level closed flow 
corresponds to choosing one of these blocks; the relevant colour algebra at the real-emission 
level is then that of the block so determined. In other words, when performing an MC sum 
over Born flows, each seed is associated with a block in the matrix of real colour flows; 
within this block, the colour algebra is performed exactly, i.e. without using MC methods. 

Let me finally discuss the case of fixed real flows. This scheme has the nice property 
that the colour-flow matrix can be factored out of the subtraction procedure. In other 
words, all of the four terms in the linear combination of eq. (1.4) are proportional either to 
C(E', S) (for gluon amplitudes) or to C(T',T) (for quark-gluon amplitudes). This implies 
that colour never enters into the definition of subtraction terms if not in a trivial way, which 
does not happen (owing to the definition of colour-linked Born's) in the other formulations 
of the FKS subtraction given here, let alone in the colour-summed one. A pleasant implica- 
tion of this fact is the disappearance of the colour-linked Born's from the soft subtraction 
terms (see eqs. (2.90) and (3.105)). This, which is nothing but the colourless nature of 
the subtraction procedure, together with the possibility of defining the soft and collinear 
kinematics so as the underlying Born kinematics coincide, implies the following equation 
(which I write for quark-gluon amplitudes, for greater generality, and without considering 
the case of the g — >• qq branching, which I shall discuss later): 



where the arrow denotes the soft, collinear, or soft-collinear limit. The kernel /Qj is thus 
able to describe the singular behaviour of the interference between dual amplitudes in these 
three limits. The only contribution to the divergences that cannot be included in /Qj is 
the azimuthal-dependent part of the collinear limit, for which the kernels /C^?t need be 




(6.1) 
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introduced, owing to the different structure of the factorized Born-level dual amplitudes 
(the helicity of the branching parton is not summed over in the last two terms on the 
r.h.s. of eq. (6.1)) 6 . As is known, the azimuthal terms vanish upon integration, but must 
nevertheless be included in the subtraction, in order to have a local cancellation of diver- 
gences in all phase-space points, and not only at the integrated level. Needless to say, in 
eq. (6.1) and in what follows all terms of collinear origin are present only if parton j is not 
a massive quark. The explicit form of the kernels that appear in eq. (6.1) can be worked 
out from eqs. (3.105) and (3.109). One obtains: 
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(6.2) 
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(6.4) 
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Several comments are in order here. First of all, being clear that the results obtained in 
sect. 3 are valid for any partons, I have relabeled n + 1 and s as i and j respectively, 
consistently with the role of Sij in eq. (6.1); it is however understood that parton i is 
a gluon. The SOFT and COLL tags in eqs. (6.2)-(6.4) imply that the four momenta 
used to compute the corresponding quantities are those relevant to the soft and collinear 
configurations respectively. In the antenna or dipole formulations, the analogues of these 
configurations are explicitly defined by means of invariants. In FKS, one does not need 
an explicit definition; any maps induced by the real-emission phase-space parametrization 
will do, provided that the latter satisfies the condition that the underlying Born kinematics 
obtained by taking the soft and collinear limits coincide 7 . The exactly collinear configura- 
tions used in the computation of the second term on the r.h.s. of eq. (6.2), and in eqs. (6.3) 
and (6.4), is the reason why no S function appear there - in such a limit, Sij = 1 [13]. 
The (+) index attached to the Altarelli-Parisi kernel in the second term on the r.h.s. of 
eq. (6.2) implies that such a kernel must be understood as being subtracted by means of 
a plus prescription. It is easy to convince oneself (see e.g. eq. (2.29)) that this is sufficient 
to take into account the soft-collinear limit (last term in eq. (1.4)). Although in principle 
the same should be done for the azimuthal kernels of eqs. (6.3) and (6.4), in practice this 
is not necessary, owing to the fact that Q a b*( z ) vanish in the soft limit. I finally point out 
that eqs. (6.1)-(6.4) are valid also in the case of g — > qq branching, with only changes in 



Alternatively, one can define in a straightforward manner a unique kernel, as a tensor in the space of 
the helicities of the branching parton. This is what is done e.g. in dipole subtraction. 

7 As was mentioned at the end of sect. 2.1, this condition is necessary for eq. (6.1) to hold. It is however 
not mandatory in general for FKS subtraction, but has been found to be convenient from the numerical 
point of view. 
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notation. However, because of the vanishing of the soft limit in this case, the final result 
would be identical to that already presented in eqs. (3.156) and (3.162)-(3.165). 

I have already shown in sect. 5 that Born- like contributions to the NLO cross section 
can also be cast in the form of terms with fixed real flows. As far as the Born proper is 
concerned, in order to integrate it simultaneously with the other terms one can use the 
same trick as that used to transform eq. (5.2) into eq. (5.8) - this is possible since the Born 
can be split into several components associated with Sij functions, as explained in ref. [11], 
thus effectively defining Born contributions at fixed real flows. Another possibility is that 
of exploiting the block structure of the real colour-flow matrix discussed above: for a given 
real flow T, there exists a Born flow 7 such that r € £(7). In this way, for any random 
choice of T in the MC sum over real flows, one also performs an MC sum over Born flows 
(which is unbiased, at least as long as T is chosen flat). 

I have presented the results in the form of unphysical — > m processes for ease of 
notation. It is clear that they have the same structure, colour-wise, when two of the final- 
state particles are crossed into the initial state. It should be kept in mind that for physical 
2 — > m — 2 processes there is contribution to the NLO cross section which is essentially a 
Born-like one, except for the fact that it features an additional integration variable. This 
contribution, called "degenerate (n+l)-body" in ref. [11], has the same form as the collinear 
limits discussed at length in this paper, and can therefore be manipulated in exactly the 
same way. 

Antenna and dipole subtractions. It seems appropriate to conclude this discussion by 
stressing the connections between FKS, antenna, and dipole subtractions, which are best 
uncovered when a fixed-real-flow scheme is adopted in the former. Firstly, given the fact 
that antenna subtraction is formulated in terms of dual amplitudes or of quantities closely 
related to those, it is natural that FKS at fixed real flows use the same matrix elements 
limits as in antenna (see e.g. eq. (3.105) and its counterpart obtained in ref. [32], in the 
context of a process-specific computation). As was mentioned before, this is reassuring, 
since the factorization formulae one starts from do not coincide in the two approaches. 
From the kinematical point of view, however, FKS and antenna are still fairly different. 
In particular, in FKS one of the two radiators that define an antenna is not needed (the 
other radiator may be identified with the FKS sister), and this is reflected in the different 
choices made for the phase-space parametrizations (which in FKS do not depend on colour 
connections) . 

Secondly, the methods used here to manipulate the colour-linked Born's (which for 
example lead one from eqs. (3.85)-(3.87) to eqs. (3.104)-(3.106)) can also be used to carry 
out the same operations on dipoles. More specifically, the structure of the colour kernel of 
a dipole X>ij,fc is (in the notation of this paper) Q b (ij)Q b (k), up to an overall Casimir. 
The most involved case is when either i or j (or both) is a gluon - be it i just to fix 
the notation. Then, one can identify ij and k of the colour kernel above with I and k 
in eq. (3.86), i with the FKS parton, and proceed as done in this paper 8 . In this way, a 

8 When j is a gluon, T>ij,k also subtracts the soft singularities associated with j, whereas in FKS the soft 
singularities of the FKS sister are damped. In order to simplify the discussion, I assume here that such a 
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linear combination of dipoles will emerge, analogue to that featuring the eikonal factors in 
eq. (6.2): 

£ E a P v a* {i r/ ( r ' _1 « ->• rcr- 1 ^) e/?)) . (6.5) 

a=-l,l/3=-l,l 

As was the case for the colour-linked Born's, all the dipoles in eq. (6.5) will factorize the 
same colour-flow matrix element C(T',T). However, different dipoles will factorize dual 
amplitudes (which sandwich helicity-dependent kernels) computed with different kinemat- 
ics, and therefore the analogue of eq. (6.1) cannot be written in this case. This problem can 
be avoided by using variants of the dipole formalism that curb the proliferation of reduced 
kinematics (see e.g. refs [33, 34]). 

The similarities among the various subtraction formalisms need not be surprising; as 
was shown in ref. [35], all of them must have the same underlying structure, and differences 
arise when choices are made for the projections that map resolved kinematic configurations 
onto unresolved ones, and for the definitions of subtraction terms away from the zero- 
measure soft and collinear regions. Colours somehow blur the picture, which becomes clear 
again if one works in schemes such as the fixed-real-flow one (which may not be unique in 
this respect). 

7. Conclusions 

The results of this paper will allow the implementation of the FKS formalism in a colour- 
friendly way, since the elementary ingredients of the subtraction procedure can be defined 
at fixed colour configurations or colour flows. This is a necessary condition in order to 
be able to perform the colour algebra with Monte Carlo methods, and thus to tackle the 
computation of large-multiplicity processes. The formulation of the subtraction at fixed 
colour flows can be used to organize the calculations as a systematic expansion in 1/N, 
since all results given here are exact to all orders in N. When the flows are fixed at the 
real-emission level, the colour and Lorentz structures completely decouple, and this allows 
one to define colourless kernels that can simultaneously describe the soft, collinear, or soft- 
collinear behaviour of the matrix element squared. In this scheme, one sees more clearly 
the connections between the FKS subtraction method and the antenna and dipoles ones, 
which I have briefly discussed. 
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damping is performed also in dipoles, thus a-symmetrizing the roles of i and j. Many other options are 
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A. Conventions for colour matrices 



a,b=l 
N 2 — l 



Contrary to the usual conventions, in this paper I normalize the Gell-Mann matrices in the 
same way as the SU(iV) generators. Hence: 

N 2 -l 

[X a ,X b ] = i £ f abc X c , (A.l) 

c=l 

Tr (A a A 6 ) = T F 5 ab , (A.2) 

N 2 -l 

{X a \ a ) t3 = C F 5 l3 , (A.3) 

a=l 
N 2 -l 

£ f abc f abd = C A 5 cd , (A.4) 
=1 

v 2 — 1 

a=l ^ ' 

with the usual colour factors 

T F = \, (A.6) 
AT 2 - 1 

C ' = ^r- (A - 7) 

C A = A^. (A.8) 
B. Colour operators, flows, and colour conservation 

As shown in sects. 2.1 and 3.1 for gluon and quark-gluon amplitudes respectively, the 
consistency between the soft and collinear limits of the matrix elements follows from the 
colour-conservation identities of eqs. (2.45) and (3.47). These identities are in turn a 
consequence of expressing the scattering amplitudes in terms of dual amplitudes, as done 
in eqs. (2.4) and (3.14), as I shall show in this appendix. 

An important by-product of this proof is the fact that the colour operators Q b , while 
defined in a natural way in the spaces of colour configurations as in eqs. (2.33)-(2.34) 
and (3.36)-(3.38), can also be easily interpreted in terms of flows. I shall therefore introduce 
operators acting on flows that have an explicit correspondence with Q b , and which will be 
instrumental in deriving a formulation of the FKS subtraction alternative to that at fixed 
colour configurations. 

B.l Gluon amplitudes 

Since eq. (2.45) is an identity in the colour space of n + 1 gluons, it can be rewritten as 
follows 

n 

X = <oi, • • ■ a n+1 \Q b (k)\A^) = V {ai}^ 1 . (B.l) 
k=l 
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Note that eq. (B.l) is equivalent to saying that colour conservation applies to (n + l)-gluon 
amplitudes at fixed colour configurations. Using eq. (2.13), one can rewrite 

n 

X a = J2(^--^n+i\Q b (k)\A^(a)). (B.3) 



k=i 



Therefore, if one can prove that X a = for an arbitrary a, eq. (B.l) will follow (the 
converse is obviously not true, and thus X a = is a stronger conditions than X = 0). By 
using eqs. (2.11), (2.33), and (2.34), one obtains 

X a = Y a A^(a) , (B.4) 

where 

n 

Y a = ^2^2(T a " +1 )a k b k Tr^A"^ 1 ' . . . A aCT ( CT ~ 1 ( fc )- 1 )A 6 ' i: A aCT ( CT ~ 1 «+ 1 ) . . . \ a °( n ^ (B.5) 

k=l b k 
n 

= X^ Tr (^ a<7(1) • • • A" 1 ^ 1 ^)- 1 ) A afe ,A a ™ +1 \ a "(*-Hk)+i) . . . \ a *(n)^ . (B.6) 
k=i 

Hence, in order to prove that X a = 0, one must prove that Y a = 0. I start from observing 
that eq. (B.6) is a linear combination of traces obtained by inserting A a ™ +1 into the traces 
that appear in |^.( n )(cr)). Such insertions can be conveniently represented in terms of flows. 
This can be done by defining the following operators: 

I+(i)a = I + (i)(a(l),...a(n)) = (a(l), . . . a(i), n + 1, . . . a{n)) , (B.7) 
J_(i)<7 = /-(i)(<7(l), . . . a(n)) = (<7(1), . . . n + 1, <r(i), . . . a(n)) . (B.8) 

In other words, inserts the number n + 1 into the list defined by a after 

(before) the i th member of the list. One can now rewrite eq. (B.6) as follows, using the 
shorthand notation of eq. (2.5): 

n 

Y a = £ {A (l + (a-Hk))a) - A (/-(^WV) } (B.9) 
k=i 

A comparison of eq. (B.9) with eq. (B.3) shows the relationship between the Q b and the 
I± operators: 

Q b <— >• /+-/_. (B.10) 

In keeping with the fact that Q b creates a colour state associated with gluon n + 1, the 
operators I± transform a n-gluon flow into an (n + l)-gluon flow. A technical difference 
between the operators Q b and I± is that the argument of the former is a particle label, 
whereas the arguments of the latter are the positions in the list of particles that defines 
the flow. This is justified by the fact that when manipulating flows it is more convenient 
to deal with positions rather than particle labels - singularity configurations are due to 
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Furthermore 



particles that are adjacent in a flow. It is obvious that summing over particle labels is 
equivalent to summing over positions in a flow. Hence, from eq. (B.9): 

n 

Y a = {A (I+(i)a) - A (J_(i)<7) } (B.ll) 

n n 

5>(J_(i)<7) = A(/_(lV) + E A ( 7 -(^) 

i=l i=2 

n 

= A(J + (n)<7) + 2>(J + (i-l)<r) 

i=2 

= A(/ + (nV) + ^A(/ + (iV) 
i=i 

n 

= 2>(J + (i)<7), (B.12) 
i=i 

where the various manipulations follow from the invariance of the trace under cyclic per- 
mutations, and from the property 

J_(i)<7 = I+(i-l)<7, 2<i<n, (B.13) 

which trivially follows from the definitions of I±(i)- By replacing eq. (B.12) into eq. (B.ll) 
one proves that Y a = 0. As discussed at the beginning of this section, this not only proves 
eq. (2.45), but also that 

n 

E#)M (n) W>=0' (B-14) 
k=i 

In other words, colour conservation also holds at fixed Born flows; this is far from surprising 
from the physics viewpoint. 

B.2 Quark-gluon amplitudes 

In the case of quark-gluon amplitudes, the proof of eq. (3.47) proceeds through proving 
the analogue of eq. (B.l), i.e.: 

n 

X= (a-2 q ,...a n+1 \Q b (k)\A^) = V {a^ 2? . (B.15) 

k=-2q 

Following what is done in the case of gluon amplitudes, one writes 

X= £ X,, (B.16) 

n 

X,= £ (a- 2q ,...a n+1 \Q b (k)\A^ n \^). (B.17) 

k=-2q 
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As was done before, I shall show that the condition X 7 = holds, thereby proving eq. (B.15) 
as well. Using eq. (3.29), one gets 



x 7 = y 7 ^)( 7 ), 



(B.18) 



where now 



y 7 = Yl A ({cJ,7) Yl (a-2 q ,...a n+1 \Q b (k)\c„ 2q ,...c n ). (B.19) 

Owing to the factorized form of the colour structure, eq. (3.15), one can associate a colour 
structure with each colour line: 

A ({a,}, 7 ) = N~M f[A({a t }, lp ) , 
P =i 



A({aJ, 7p ) = (A aCT( *p-i +1) ...A a ^>) 

V / a— v a 



(B.20) 
(B.21) 



P a il( — p-q) 



Using the definition of the operator Q b (k), eq. (B.19) then becomes: 



Y 7 = AT-p(t) Y 



p =i 



II A 



7r 



. r=l 



\ 



y 



(B.22) 
(B.23) 



Equation (B.23) states formally that the action of the operators Q h onto the scattering 
amplitude can be conveniently rewritten colour line per colour line. It also implies that in 
order to prove I 7 = 0a sufficient condition is that of proving Y J:P = for all p, which is 
what I set out to do now. First, by using eqs. (2.34), (3.37) and (3.38), one obtains: 

E (® an+1 A (WiWk, c k , 7 p) = 



f\ a ^( t p-i+ 1 ) . . . ^V'W-i) A afc ,A° n+1 A aCT ( CT_1 ( fc )+ 1 ) . . . A^Cf') , (B.24) 

EW a " +I WU A (KU,cM P ) = (a^A-^-i+D ... AM«p)) , (B.25) 



P a M(-f>-9) 

E^ a " +1 (^))a fcCfc A({a l }^ fc , Cfc , 7p ) = -(A a ^-^» ... A a *W"+i) , (B.26) 

for the cases when A; is a gluon, a quark, or an antiquark respectively. Note that in 
the latter two cases the requirement that k 6 7 P implies A; = — p and fc = n(—p — q) 
respectively. Precisely as in the case of gluon amplitudes, eqs. (B.24)-(B.26) suggest the 
use of operators acting on flows. For consistency with what was done before, the arguments 
of these operators will have to be the positions of the particles that appear in flows. My 
conventions for such positions are the following. The colour line of eq. (3.26) 



7 P 



-P ; <r(t p -i + 1), . . . cr(t p ); fi(-p - qfj 



(B.27) 
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corresponds to positions: 



(-p;t p _i + l,...t p ;-p-g). 



(B.28) 



This implies that — p is adjacent to t p _i + 1, and that —p—q is adjacent to t p . It is therefore 
convenient to define the operations © 1 and 1, that will serve to move across the list in 
eq. (B.28): 



i © 1 

i © 



1 = < 

1 
-1 



tp-i + 1 
i + 1 
-p-q 

-p 
i-l 

i G 1 © 1 = i. 
iel. 



i = -p, 

tp-i + 1 < i < t p , 
% — tp , 

i = tp-i + 1 , 
tp-i + 1 < i < t p , 
i = -p-q, 



(B.29) 



(B.30) 

(B.31) 
(B.32) 



The definitions in eqs. (B.29) and (B.30) imply that one can move continuously across a 
given colour line, but cannot pass continuously from line j p to line j p ±i- This is consis- 
tent with the physical interpretation of colour lines, which from the colour viewpoint are 
disconnected from each other, and with the fact that when studying singularities at fixed 
flows the emphasis is on adjacent particles. Equations (B.28)-(B.30) are meant to hold for 
all colour lines 7 P belonging to a given flow 7. Note that in the case of eq. (3.23), i.e. when 
there are no gluons on colour line 7 p , — p and — p — q are contiguous, and therefore 



(-p) el=-p-q 

(-p-q) ei = -p 



if tp— 1 — tp . 



(B.33) 



Finally, I shall denote by 7~ 1 (/c) the position of particle k in flow 7 (exactly as cr _1 (A;) 
denotes the position of gluon k in flow a). Thus, eqs. (B.24)-(B.26) can be rewritten as 
follows: 



Cfe 



A({a t },I + ( 7 -\k)) 7p ) - A({a,},/-(7~ 1 «)7 P ) , 
E^ an+1 «)a fcCfc A({a,}^, Cfc ,7 P ) = A({a,},/ + (7- 1 (fc))7p) 



Cfc 



E(Q an+1 «W fc A c fc)7p ) = -A({a i },/47" 1 W)7p) 



(B.34) 
(B.35) 

(B.36) 



Cfe 



As in the case of gluon amplitudes, the result of acting on the list associated 

with the flow is that of inserting the number n + 1 after (before) the i th member of the list. 
Given the factorized form of the flow, eq. (3.25), the operators I±(i) can be equivalently 
understood as acting on the colour flow 7, or on the colour line 7 p which includes the i th 
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position. The conventions adopted before imply that: 



I + (r\k)) =I+(-p) k = -p, (B.37) 

/-(7 -1 (*0) = I-(-p ~q) k = fi(-p - q) . (B.38) 

Putting all this together, one obtains: 

tp tp © 1 

Y 7 , p = £ A({aJ,/ + (i) 7p ) - J2 A({aJ,/_(i) 7p ) . (B.39) 

i=t p -i+101 i=t p _i+l 

By construction, the analogue of eq. (B.13) holds: 

i-(i)7 P = /+(«el)7 P , tp-i + 1 < i < t p 1. (B.40) 

Hence 

tp © 1 tp©l 

E A ({<*}, J_(i) 7p ) = E A({a 4 },/ + (;el) 7p ) (B.41) 

i=t p _ 1 +l Z=t p _i + 1 

tp 

= £ A ({a,}, /+(*)>) , (B.42) 
i=t p -i+iei 

where in the last equation I relabeled the sum variable i — >■ i0l. By plugging eq. (B.42) into 
eq. (B.39) one gets finally gets Y 7iP = 0. It is easy to convince oneself that the derivation 
above also applies to the case of a colour line with no gluons attached, since in such a case: 

I+(-p)lp = I-(~P ~ q)l P if tp-i = t p . (B.43) 
To take this fact into account, one may perform the following formal replacement: 

tp tp 

E /« E (B.44) 

i=t p -i+iei i=t p _i+iei 

In practice, the replacement in eq. (B.44) will always be understood. Given that the 
arguments above apply to any p, one has indeed proved eq. (3.47), and also that 

E Q b (k)\A^^)) = 0, (B.45) 

k=-2q 

analogously to what happens in the case of gluon amplitudes. 
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